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Abstract 

Plant growth is continuous and modular, a combination that allows morphogenesis by cell division and elongation and 
serves to facilitate adaptation to changing environments. The pleiotropic phenotypes of the harlequin (hlq) mutant, 
isolated on the basis of ectopic expression of the abscisic acid (ABA)- and auxin-inducible proDc3:GUS reporter gene, 
were previously characterized. Mutants are skotomorphogenic, have deformed and collapsed epidermal cells which 
accumulate callose and starch, cell walls abundant in pectins and cell wall proteins, and abnormal and reduced root 
hairs and leaf trichomes. hlq and two additional alleles that vary in their phenotypic severity of starch accumulation 
in the light and dark have been isolated, and it is shown that they are alleles of b//i3/ftyp6/rh/3/Topoisomerase6B. 
Mutants and inhibitors affecting the cell wall phenocopy several of the traits displayed in hlq. A microarray analysis 
was performed, and coordinated expression of physically adjacent pairs/sets of genes was observed in hlq, suggest- 
ing a direct effect on chromatin. Histones, WRKY and IAA/AUX transcription factors, aquaporins, and components 
of ubiquitin-E3-ligase-mediated proteolysis, and ABA or biotic stress response markers as well as proteins involved 
in cellular processes affecting carbon partitioning into secondary metabolites were also identified. A comparative 
analysis was performed of the hlq transcriptome with other previously published TopoVI mutant transcriptomes, 
namely bin3, bin5, and caa39 mutants, and limited concordance between data sets was found, suggesting indirect 
or genotype-specific effects. The results shed light on the molecular mechanisms underlying the det/cop/fus-\\ke 
pleiotropic phenotypes of hlq and support a broader role for TopoVI regulation of chromatin remodelling to mediate 
development in response to environmental and hormonal signals. 
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Introduction 

Plants are sessile and have evolved an open pattern of growth 
and development to give plasticity to their responses to 
changing environments. A central paradox of plant biol- 
ogy is that there exist only a handful of compounds (phy- 
tohormones) that mediate many complex processes in plant 
growth. A related conundrum is that, despite the elucida- 
tion of hormone structures, their biosynthetic pathways, 
and cognate hormone receptors, an understanding of the 
mechanisms of hormone action and 'sensitivity' remains elu- 
sive. Abscisic acid (ABA) is a plant hormone that mediates 
a myriad of physiological processes in growth and develop- 
ment, including cell division, water use efficiency, and gene 
expression during seed development and in response to envi- 
ronmental stresses, such as drought, chilling, salt, pathogen 
attack, UV, and high light (Rock et al, 2010). It was previ- 
ously shown that the uidA (^-glucuronidase; GUS) reporter, 
under the transcriptional regulation of the carrot (Daucus 
carota) Late Embryogenesis- Abundant Dc3 promoter in 
transgenic Arabidopsis thaliana seedlings, is ABA and auxin 
inducible within the root zone of elongation and within the 
vasculature, and inducible in guard cells by ABA (Chak et al. , 
2000; Rock and Sun, 2005). The abil-1 and abi2-l mutations 
reduce ABA-inducible proDc3:GUS expression in root tis- 
sues and interact with (attenuate) auxin signalling leading 
to proDc3:GUS expression. Furthermore, the auxl and axr4 
mutants, which are affected in auxin transport/homeostasis, 
show a hypomorphic effect on ABA-inducible proDc3:GUS 
expression, demonstrating that ABA and IAA signalling 
pathways interact in roots (Rock and Sun, 2005). 

Several reports have established links between ABA and 
the cell wall. For example, ABA-dependent callose accumula- 
tion is required for |3-aminobutyric acid-induced resistance 
to the pathogens Altemaria brassicola and Plectosphaerella 
cumumerina (Ton and Mauch-Mani, 2004; Ton et al, 2005). 
ABSCISIC ACID (ABA)-INSENSITIVE-8 is allelic to 
ELONGATION DEFECTIVE- 1 and KOBITO-1, which 
establishes a link between ABA and the cell wall because muta- 
tions at this locus disrupt ABA-regulated gene expression, 
sugar sensitivity, cell elongation, cellulose synthesis, vascular 
differentiation, and root meristem maintenance (Cheng et al , 
2000; Pagant et al , 2002; Brocard-Gifford et al , 2004). ABI8/ 
KOB1/ELD has recently been shown to regulate movement 
of transcription factors through plasmodesmata to modu- 
late stomatal patterning (Kong et al, 2012). Mutations in the 
cellulose synthase gene (CESA8)I IRREGULAR XYLEM1 
result in resistance against drought, salt, and osmotic stresses, 
and in elevated ABA-inducible gene expression (Chen et al , 
2005). 

Because there is no cell migration in plants, and cell walls 
are formed concomitant with cell division, the timing and ori- 
entation of cell division and expansion, mediated ultimately 



by hormone signalling to the cell wall, are the primary forces 
that shape plants. Identification and cloning of mutants that 
affect plant shape and environmental responses can shed 
light on the molecular mechanisms that control morpho- 
genesis and hormone action in response to stresses. So far, 
it has proved difficult to predict the nature of the gene prod- 
ucts involved in hormone sensitivity. Previously the harlequin 
{hlq) mutant was isolated based on ectopic expression of the 
proDc3:GUS transgene (Subramanian et al, 2002). Here it is 
reported that hlq, which is de-etiolated in the dark with open 
cotyledons and short hypocotyls, encodes an allele of brassi- 
nosteroid insensitive3lhypocotyl6lroot hairless3 (bin3lhyp6l 
rhlJ), the plant-specific homologue of type lib archaebacte- 
rial topoisomerase VI subunit B. 

Type II topoisomerases have a catalytic cycle producing dou- 
ble-stranded breaks (as distinct from type I single-stranded 
breaks) that are associated with the occurrence of transient 
DNA-topoisomerase covalent complexes implicated in DNA 
replication, recombination, chromosome segregation, chro- 
matin remodelling, and transcription (Varga-Weisz et al, 
1997; Vos et al , 201 1). Type II topoisomerases can act as tran- 
scriptional repressors, by binding to promoters and blocking 
the formation of stable pre-initiation complexes, which can be 
relieved by the addition of sequence-specific transcriptional 
activators (Brou et al. ,1993; McNamara et al. , 2008). Topo6B 
has conserved ATP binding and hydrolysis domains, and a 
conserved basic motif (B4) of unknown function (Hartung 
and Puchta, 2001). In addition to the archaeal TOP6A/BIN5/ 
RHL2/SPOH-3 and TOP6B subunits, Arabidopsis TopoVI 
function requires the activity (Simkova et al, 2012) of two 
small subunits, Atlg48380 IRHL1 (Schneider et al, 1998)/ 
HYP7 (Sugimoto-Shirasu et al, 2005) and At5g24630IBIN4 
(Breuer et al, 2007)/ MIDGET (MID) (Kirik et al, 2007). 
The Arabidopsis BIN4/MID subunit of TopoVI can interact 
in a yeast two-hybrid assay with TFIIB, which is involved in 
RNA polymerase II recruitment and transcription initiation 
in eukaryotes (Evans-Roberts et al, 2010; Szklarczyk et al, 
2011). Molecular characterization of mutants of bin3lhyp6l 
rhl3ltop6b and the regulatory subunits rhlllhyp7 and bin4l 
mid that bind each other, to TOP6A, and to COP1 (Schrader 
et al, 2013) led to claims of TopoVI as having complex roles 
in growth and development. Namely, it is essentially involved 
in brassinosteroid sensitivity (Yin et al, 2002) but probably 
not hormone cross-talk (Nemhauser et al, 2006), in decat- 
enation of endoreduplicated chromosomes throughout the 
nucleus (Sugimoto-Shirasu et al., 2002), in DNA replica- 
tion/repair and coupling to cell cycle arrest that controls cell 
expansion and proliferation (Hartung et al., 2002; Schrader 
et al, 2013), including in response to nematode parasites 
(Vieira et al, 2013). Moreover, TopoVI is involved in tran- 
scriptional silencing through chromatin organization (Kirik 
et al, 2007), and as a genuine component of singlet oxygen 
retrograde signalling to the nucleus, being both a positive 
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and a negative integrator of different reactive oxygen species 
(ROS) signals and environmental cues (Simkova et al, 2012). 
Further studies have shown that the constitutive expression 
of rice OsTOP6A or OsTOP6B increases the expression of 
stress-responsive genes, and confers abiotic stress tolerance 
to transgenic Arabidopsis plants (Jain et al., 2006, 2008). 
The present phenotypic characterization and transcriptome 
profiling results for hlq seedlings now demonstrate a funda- 
mental role for TopoVI in plant gene regulation, and provide 
additional insight into the molecular mechanisms and pro- 
cesses of TopoVI which integrate environmental (e.g. light 
and dark) and internal (e.g. hormone) signals via chromatin 
remodelling to control plant growth and development. 

Materials and methods 

Plant materials 

Col-0 (stock # CS60000), abi2-l (CS23), Ler (CS20), prcl-1 
(CS297), botl-l (Bichet et al, 2001), korl-1 (CS298), and T-DNA 
insertion lines SALK_024455C (subsequently referred to as the 
hlq-3 allele) and SALK_1 40704 (hlq-2) were obtained from the 
Arabidopsis Biological Resource Center (Ohio State University, 
Columbus OH, USA; http://abrc.osu.edu/). The hlq-ll+ (CS68702), 
F 2 mapping stocks (CS68703, CS68704), hlq-2l+ (CS68705), and 
hlq-3/+ (CS68706) seed stocks have been deposited at the ABRC. 
Seeds were sown on plates containing 0.5x Murashige and Skoog 
salts (Research Products International, Mt. Prospect, IL, USA), 
0.8% sucrose, and 0.5% phytagel (Sigma-Aldrich, St. Louis, MO, 
USA). Plates were kept in the dark at 4 °C for 3 d before being trans- 
ferred into a growth chamber. The growth conditions were 21 °C 
with continuous light (~ 100 uF nr 2 s~'), except for starch quantita- 
tion assays where plants were given a diurnal light cycle of 12 h light 
and 12 h dark, and roots and shoots (four biological replicates of 
~10mg each) were harvested at the end of the respective cycle. Plates 
were kept in the chamber for 5 d before scoring for the hlq mutant 
morphology. 

Genetic mapping 

Previously, 5100 Mj cibi2labi2 homozygous plants of a line that 
carried two independent proDc3.GUS reporter genes were muta- 
genized with ethylmethane sulphonate (EMS), and M 2 clonal lines 
were screened for ABA- or auxin-inducible GUS expression in roots. 
These studies resulted in isolation of a new allele of rootylsuperrootll 
hookless3 (Sun, 2003) and two single gene nuclear mutants: hlq and 
short blue root (sbr) (Subramanian et al, 2002). These mutants have 
novel ABA- and auxin-inducible proDc3:GUS gene expression phe- 
notypes attributable to ABA and indole acetic acid (IAA) responses, 
as well as pleiotropic growth phenotypes that suggest a link to 
ABA and/or auxins. The hlql+ mutant stock (subsequently termed 
hlq-1 allele) were backcrossed to the parental line abi2 three times 
to remove extraneous mutational load caused by the original EMS 
treatment and to provide supporting evidence that the pleiotropic 
phenotypes of hlq were due to a single mutation. Due to seedling 
lethality and sterility, hlql+ heterozygous plants (in the abi2-l mutant 
background) of the Landsberg erecta (Ler) ecotype were crossed with 
the Columbia (Col-0) ecotype. F 2 seeds from a single Fj plant consti- 
tuted one mapping population. Two F 2 lines that segregated for hlq 
were used for mapping experiments. The segregation in the F 2 popu- 
lation for 3: 1 wild type to hlq showed that the hlq phenotype was not 
influenced by the ecotype. No obvious differences in the hlq mutant 
morphology compared with the Ler parental background suggested 
that penetrance of the phenotype was complete. DNA was extracted 
from 2- to 3-week-old mutant seedlings as described (Xin and Chen, 
2012). The hlq recombinants were allowed to grow for 3^1 weeks to 



maximize template DNA yield. The DNA from individual recom- 
binants was used in simple sequence length polymorphism (SSLP) 
(Bell and Ecker, 1994) and cleaved amplified polymorphic sequence 
(CAPS) PCR analyses (Konieczny and Ausubel, 1993), or amplified 
using primers flanking single nucleotide polymorphisms (SNPs) iden- 
tified from the publicly available Monsanto/Cereon Arabidopsis Ler 
draft sequence (Jander et al, 2002) [available for download at The 
Arabidopsis Information Resource (TAIR) (www.arabidopsis.org)]. 
Repetitive sequences were scanned using the REPEATMASKER 
program (www.repeatmasker.org). Oligonucleotide primers were 
designed using 'Perlprimer (http://perlprimer.sourceforge.net/) 
and synthesized commercially (Sigma; http://sigma.com) and used 
for PCR amplification using PrimeSTAR HS DNA Polymerase 
(TAKARA, Shiga, Japan). Amplicons of 1 00-200 bp were targeted to 
reveal the longest repeats and hence the high probability of polymor- 
phism in agarose gels. Among the different kinds of repeats scanned, 
TA, CT/CTT, and GA/GAA repeats showed highest abundances and 
successes for developing useful polymorphic markers at 24, 22, and 
14%, respectively, of the total scanned. Information about the posi- 
tion and primer pairs for these markers is given in Supplementary 
Table S6 available at JXB online. Amplicons were desalted with a 
Qiaquick PCR Purification Kit (Qiagen, Gaithersburg, MD, USA) 
and sequenced by the Sanger dideoxynucleotide chain termination 
method by Beckman Coulter Genomics (Danvers, MA, USA; http:// 
www. beckmangenomics. com/) . 

Whole- genome resequencing of the hlq/+ heterozygote 

Genomic DNA was prepared with a Qiagen DNeasy Plant Minikit 
(Qiagen) according to the manufacturer's protocol. Paired-end 100 
base Illumina HighSeq reads as fastq files (including quality scores) 
were mapped to the Col-0 reference genome (TAIR9) using Burrows- 
Wheeler Aligner (BWA; http://bio-bwa.sourceforge.net/) software 
with default settings, except the extension gap was set to '10' (Li and 
Durbin, 2010). This resulted in a slight improvement to 45.7% of all 
reads successfully mapping to the reference Col-0 genome, and 38.7% 
of all reads properly paired, yielding an average of 25-fold coverage 
over all five chromosomes. SAMTools (http://samtools.sourceforge. 
net/) (Li et al. , 2009) was used to index and align the reads to the 
reference genome, and SnpEff (http://snpeff.sourceforge.net) was 
used to call SNPs and INDELs (Cingolani et al, 2012). There were 
223 homozygous Ler SNPs called in the — 1 12 kbp interval of inter- 
est where hlq mapped in the region chr3: 7213133. .7325482, plus 19 
deletions and 18 insertions. There were 11 heterozygous Ler SNPs 
called, one heterozygous insertion, and one heterozygous deletion. 
A published whole-genome assembly of Ler with > 320-fold coverage 
(Schneeberger et al , 201 1) was used to validate the SNP calling meth- 
odology. Both heterozygous INDELs were called correctly, as were 
92% of the Ler homozygous SNPs/INDELs in the interval (data not 
shown), whereas an additional 15 SNPs (7%) called by Schneeberger 
et al. (201 1) were not called in the present analysis. It is possible that 
the extra homozygous SNPs in the data set here and those described 
but not found in the data are polymorphic within Ler accessions. 
Notwithstanding, the high concordance between the present results 
and the published Ler SNPs validated the method. 



GUS assays 

Plant samples were developed at 37 °C for 12 h with a 1 mM solution 
of the indigogenic GUS substrate 5-bromo-4-chloro-3-indolyl (3-d- 
glucuronide (X-Gluc; Rose Scientific, Edmonton, Alberta, Canada) 
in 50mM KH 2 P0 4 (pH 7.0), 0.1 mM EDTA, 0.5mM ferricyanide, 
0.5 mM ferrocyanide, 0.05% sodium azide, 0.1% Triton X-100. After 
staining overnight, samples were immersed in 70% ethanol over- 
night. Samples were then mounted in 30% glycerol on microscope 
slides for pictures to be taken with a dissecting microscope, or were 
imaged with an Olympus BX41 light microscope for high magnifi- 
cation detailed imaging using QCapture software (v2.68.6, Silicon 
Graphics, Fremont, CA, USA). 
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Analysis of gene expression and genotypes with PCR 

Total RNA from roots or leaves of 10-day-old seedlings, or from 
whole seedlings, was extracted with Iso-RNA Lysis Reagent (5 
PRIME, Gaithersburg, MD, USA). For reverse transcription-PCR 
(RT-PCR), 2 ng of total RNA was digested with RQ1 RNase-free 
DNase (Promega, Madison, WI, USA) and reverse transcribed by 
MMLV reverse transcriptase (Promega) with Anchored Oligo-dT 
(ABgene; Surrey, UK). A cycle number of 30 was used for PCR 
experiments in Fig. 4B; for Supplementary Fig. S7 at JXB online, 28 
cycles were used for SALK_140704, and 3 1 cycles for SALK_024455. 
PCR conditions were: 55 °C annealing temperature, 1.5min exten- 
sion time, 72 °C extension temperature. PCR products were run 
on an agarose gel with a DNA 1 kbp ladder (Axygen Biosciences, 
Union City, CA, USA; www.corning.com/axygen/) for calculation 
of band sizes by linear regression. Primers used in this study are 
listed in Supplementary Table S6. 

RNA gel blot hybridization 

For RNA gel blot hybridization, samples of 20 (ig of total RNA 
were resolved on a 1.2% denaturing agarose gel and blotted to a 
Hybond-N+ membrane (GE Healthcare, Piscataway, NJ, USA) 
according to the supplier's protocol. Ambion Millenium RNA mark- 
ers (GE Healthcare) were run as controls for calculation of band sizes 
by linear regression. At3g20780 template was amplified using the 
primers listed in Supplementary Table S6 at JXB online from a cDNA 
library and was eluted after excision of the band from an agarose gel. 
Probes were synthesized using the Random Primer DNA Labeling 
Kit Ver.2 (TAKARA) with [a- 32 P]dCTP (PerkinElmer, Waltham, 
MA, USA). Hybridization was carried out with the PerfectHyb 
Plus hybridization buffer (Sigma-Aldrich) according to the manu- 
facturers instructions. A storage phosphor screen (GE Healthcare) 
was used for autoradiography and it was scanned using a Storm 860 
Phosphorlmager (GE Healthcare). Band intensities were quantified 
using the ImageQuant TL software (v2003, GE Healthcare). 

Microarray experiment 

Three individual pools of total RNA were extracted with Trizol 
reagent (Invitrogen, Carlsbad, CA, USA) from ~1000 seedlings (14 
d old) segregating hlq homozygotes from self-fertilization of the 
abi2-llabi2-l-hlql+ genotype after germination on Petri plates as 
described and kept vertically under continuous light at 21 °C. The 
control (wild-type) samples were three pools of normal individuals 
from the same plates and comprised 'wild type' (i.e. abi2-llabi2-l 
homozygous) and abi2-llabi2-l-hlql+ heterozygous genotypes used 
for microarray analysis. In vitro transcription was carried out using 
an Amino Allyl MessageAmp aRNA Amplification Kit (Ambion; 
cat#: AM1753) according to the manufacturer's protocol (Grand 
Island, NY, USA). A lOOOng aliquot of Cy3- or Cy5-labelled cRNA 
from different replicate genotype samples was purified separately, 
combined in equimolar amounts including dye swaps across geno- 
types, and mixed with hybridization buffer before being applied to 
three microarrays. Each sample was run on a Bioanalyzer (Agilent 
Technologies, Santa Clara, CA, USA) before and after labelling. The 
ratio of red (Cy5) and green (Cy3) fluorescence intensities for each 
spot is indicative of the relative abundance of the corresponding tar- 
get molecule in the two mRNA target samples. A fourth microarray 
was hybridized with balanced pooled samples from the same dye- 
labelled biological replicates as a technical replicate. Hybridizations 
were performed on the Operon Arabidopsis Long Oligonucleotide 
Microarrays V3 GEO platform GPL4570. Conditions of hybridi- 
zation and washing were according to the protocol published by 
Galbraith et al. (2011). After washing, microarrays were scanned 
using the GenePix Autoloader 4200AL Scanner (Molecular Devices, 
Sunnyvale, CA, USA) with laser excitation at 532 nm and 635 nm at 
the resolution of 10 pixels urn -1 , and saved as 16-bit grey scale multi- 
image TIFF files. Intensity values were extracted using GenePix 
Pro 6.0 (Molecular Devices) and saved as 'gpr' and 'txt' files for 



each block individually. The data for each array were lowess nor- 
malized followed by analysis of variance (ANOVA) of differential 
gene expression using empirical Bayes methods to moderate the 
standard deviations between genes. This was done by the software 
package Linear Models for Microarray Data (limma) (Smyth et al. , 
2005) written in R language under the BioConductor platform 
(Gentleman et al, 2004) version R2.5.1 (http://www.r-project.org/). 
Normalization 'within' and 'between arrays' scaled the log2-ratios 
of each two-colour experiment to have the same median-absolute 
deviation across arrays. A dye effect coefficient was calculated in 
addition to the genotype coefficient and signal amplitude across all 
probes and samples. Original MIAME-compliant data are stored 
at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/ 
geo/) with the locator GSE45806. The key reporters referenced in 
the literature for retrograde ROS signalling [with fold change (FC) 
observed in hlq seedlings following in parentheses] are: At5g40010l 
AAA-ATPasel (FC=1.08), At3g61190IBONSAI ASSOCIATION 
PROTEIN1 (FC= 1.68, PO.0005), and At5g01600IFERRITINl 
(FC= 1.48, P<0.03). The lack of significant up-regulation for AAA- 
ATPasel in hlq mutants is consistent with previous results for 5-day- 
old seedlings of top6alcaa39 (Simkova et al, 2012). 



Lignin, callose, pectin, and starch staining/quantitation, 
dichlobenil and tunicamycin treatments 

For staining for lignin with phloroglucinol, the seedlings were treated 
with a mixture of methanokacetic acid (9:1) to remove chlorophyll. 
The seedlings were then cleared in chloral hydrate:water (8:2) for 
2-5 min and stained with phloroglucinol (1% phloroglucinol in 
6 N HC1) for 5 min (Zhong et al, 2000). Stained seedlings were 
observed and photographed using an Axio-vert inverted microscope 
or Axiophot microscope. Inverse phase contrast was used to bet- 
ter visualize lignin staining. For staining lignin with basic fuchsin, 
the seedlings were treated with a mixture of methanokacetic acid 
(9:1) to remove chlorophyll and then cleared in 10% NaOH at 60 °C 
for 12h. Cleared seedlings were stained with 0.01% Basic Fuchsin, 
mounted in 50% glycerol, observed by epifluorescence using a pro- 
pidium iodide (PI) filter, and photographed. 

Callose staining was performed with slight modifications from 
Buogourd et al. (2000). The seedlings were treated with a mixture 
of methanokacetic acid (9:1) to remove chlorophyll and then re- 
hydrated in water. The seedlings were stained with a 0.025% aque- 
ous solution (1 in 20 dilution of 0.5% stock solution) of aniline blue 
(Sigma-Aldrich) for 30 min. The seedlings were then washed twice 
with water for 15 min each to remove excess aniline blue, mounted 
in 50% glycerol, observed by epifluorescence, and photographed 
using a 510nm filter. Callose was quantified from seedlings of differ- 
ent ages according to Kohler et al. (2000). Briefly, the seedlings were 
weighed, ethanol dehydrated and decolorized for 2 d, homogenized, 
and extracted with the loading mixture (0.1% aniline blue, 1M gly- 
cine/NaOH, 1 M HC1). The fluorescence was measured in a fluores- 
cence spectrophotometer at 393 nm excitation and 479 nm emission 
wavelengths. 

Starch staining was with Lugol's elemental iodine:potassium 
iodide (IKI; 5.7 mM iodine, 43.4 mM potassium iodide) after pig- 
ment removal and dissolution of membranes by boiling in 80% 
ethanol, as described (Kurata and Yamamoto, 1998). Starch quan- 
titation was by solubilization (boiling pulverized frozen material 
in water after extracting free glucose with 90% ethanol, 60 °C for 
5 min), complete conversion with glucan hydrolases to glucose, fol- 
lowed by enzyme-based fluorimetric assay (Abeam; www.abcam. 
com) according to the manufacturer's protocol. 

For staining pectins, the seedlings were stained with Ruthenium 
Red (SPI-Chem, Cat. No. 02603-AB) at 0.2% in water with 0.01% 
Tween for 1 h, rinsed twice in distilled water, and viewed under a 
dissection microscope. Tunicamycin or 2, 6-dichlorobenzoni- 
trile (dichlobenil, DCB; Sigma-Aldrich) was added to the growth 
medium at a concentration of 1, 2, or 3 uM. Two-day-old seedlings 
were transferred to tunicamycin-containing medium and allowed for 
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grow for a further 5 d. Their growth was monitored and phenotypes 
observed on 1 -week-old seedlings. 

Peroxidase activity 

Four-week-old wild-type and hlq mutant plants grown on agar plates 
were homogenized at 4 °C in homogenization medium [10x vol. g 
fresh weight (FW) -1 ] comprising 50 mM TRIS-HC1 (pH 7.0), 2 uM 
phenylmethylsulphonyl fluoride, and 1 |xg ml -1 of each of the pro- 
tease inhibitors leupeptin, pepstatin A, and aprotinin (Sigma). The 
homogenate was centrifuged at 3000 g for 15 min to obtain a soluble 
and crude cell wall fraction. The pellet was washed in homogeniza- 
tion medium and then incubated with 1 M KC1 (2x vol. g FW -1 ) in 
50mM TRIS-HC1 buffer (pH 7.0) for lOmin prior to centrifugation 
for 10 min at 3000 g to obtain the KC1 extract. The KC1 extract was 
microfuged at 13 000 rpm for 7 min at 4 °C to remove cell debris and 
concentrated 10-fold by ultrafitration using a 1.5ml Fugisep (lOkDa 
molecular weight cut-off membrane; Intersep). The KC1 extracts 
were equilibrated in 20mM TRIS-HC1 (pH 7.0) by diafiltration 
using Fugisep-10 membranes (Intersep) prior to loading on a 10% 
native-polyacrylamide gel (Brownleader et al , 2000). Electrophoresis 
was performed at 4 °C at 30 mA constant current. The gel was then 
immersed in a solution of 0.1% guaiacol and 0.03% H 2 0 2 in 50 mM 
potassium acetate buffer, pH 6.0 until bright blue bands appeared 
(7-10min). The reaction was stopped by immersing in 7% acetic acid. 

Microscopy 

Light microscopy was performed using a Zeiss stereomicroscope 
(Gottingen, Germany), and for micro measurements a calibrated 
microscale was used. For light microscopy of GUS-stained roots, 
the seedlings were briefly rinsed in 70% ethanol followed by sterile 
distilled water to remove excessive GUS developer and observed 
under the microscope. For 'live and dead' cell staining, seedling roots 
were embedded by allowing them to grow into 2% Phytagel minimal 
medium. Embedded roots in 1-2 mm thick slices of gel were trans- 
ferred to glass slides, immersed in water containing 100 ug ml -1 each 
of fluorescein diacetate (FDA; Sigma) and PI (Sigma), then incubated 
in the dark (15min), rinsed twice in sterile distilled water, and epifluo- 
rescence was viewed under excitation with a blue filter (450-490 nm) 
by means of a Zeiss Axiophot microscope. The sample was then 
flooded with GUS developer solution for 16h and photographed in 
bright field with the same field in view as previously documented. For 
aniline blue visualization, excitation was with UV and emission was 
visualized with a 4',6-diamidino-2-phenylindole (DAPI) filter set. 

For scanning electron microscopy, the seedlings were aligned on 
the conductive paste thinly spread over the microscope stage and 
immediately frozen under liquid nitrogen. The stage was placed 
into the cryo-chamber of the microscope and the condensed water 
removed by vacuum pumping. Then the samples were coated with 
platinum and observed under a Leica S440 scanning electron micro- 
scope (Cambridge, UK) fitted with a Fisons LT7480 cryoprep cryo- 
stage (Fisons Instruments, UK). 



Fourier transform infrared (FTIR) microspectroscopy 

FTIR has been used to 'fingerprint' the carbohydrate constituents in 
the 1200-900 cm -1 region, to detect and characterize conformational 
changes in wall components, and to determine which cross-links 
between polymers are present. Diagnostic absorbances are as follows: 
the carboxylic ester group absorbs at ~1740cm _1 , amide-stretching 
bands of protein occur at ~1650cm -1 and 1550 cm 1 , carboxylic acid 
groups on pectins absorb at 1610cm" 1 , phenolics absorb at ~1600cnr 
1 and 1500 cm -1 , and carbohydrates absorb between 1200 cm -1 and 
900cm -1 (McCann et al., 1992). IR spectra of hlq and the wild type 
were acquired with the Autoimage FT-IR micro-spectroscopy sys- 
tem (PerkinElmer; Norwalk, CT, USA). The microscope includes 
a camera and a viewing system that magnifies the visible light 
image of the sample and enables isolation of a region of interest. 



The image of the sample is displayed on the monitor visible window. 
The AutoIMAGE software (version 5.0.0 B8) enables the control of 
the operation of the microscope, and maps and collects the spectra 
from a sample. The system used a KBr beam splitter and utilized a 
liquid nitrogen-cooled MCT detector housed in the microscope. The 
Autoimage was further equipped with an automated XYZ motor- 
ized stage that was operated in the auto-focusing (Z-direction) mode 
under the software control. Both visible images and IR maps were 
obtained in the transmission mode from an area of 1 50 x 1 50 um. 
Samples of cell wall extracts (described below) were placed on a 
BaF 2 window (13 mm diameter, 2 mm thickness) and air dried at 
37 °C for 1 h. A background spectrum of clean BaF 2 surface was 
collected before each scan. The IR spectra were recorded using 128 
interferograms with 4 cm -1 spectral resolution between 4000 cm -1 and 
700 cm -1 . Twelve spectra were recorded from each sample. All FTIR 
spectra were baseline corrected and normalized. Principal compo- 
nent analysis (PCA) is a mathematical technique that is widely used 
to reduce the dimensionality of data sets. The variability in each indi- 
vidual spectrum relative to the mean of the population is represented 
as a smaller set of values (axes) termed principal components (PCs). 
This process concentrates the sources of variability in the data into 
the first few PCs, and was used in analysis of FTIR spectra from 
purified cell walls of hlq homo- and heterozygotes, and the wild type. 

Glycosyl composition analysis of cell walls 

Plant tissue (0.7-2 g FW) was collected and stored at -70 °C until pro- 
cessed. Frozen tissue was ground into fine powder in liquid N 2 with 
a mortar and pestle and extracted twice in 80% (w/w) phenokacetic 
acid:water (5:2: 1 v/v/v) to remove protein and other soluble metabolites 
(Fry, 2000). The pellet after centrifugation was washed with 70% etha- 
nol and extracted twice with 90% (v/v) dimethylsulphoxide (DMSO) 
to remove starch (Selvendran and O'Neill, 1987). Lipids were removed 
by washing the pellet twice with chloroform:methanol (2:1 v/v) and 
acetone (York et al, 1986). The cell wall fraction was dried in air and 
again in a vacuum desiccator over P 2 0 5 . The dry cell wall fraction from 
hlq mutant and wild-type plants was then analysed at the level of the 
glycosyl composition of non-cellulosic polysaccharides as previously 
described (Komalavilas et al, 1991). Briefly, cell walls were subjected 
to methanolysis in 1.5 M methanolic HC1 at 80 °C, trimethylsilylation, 
and gas-liquid chromatography (GC). For analysis including cellulosic 
polysaccharides (which are not cleaved by methanolysis), an aliquot 
of the cell wall fraction was pre-swollen in 22 N H 2 S0 4 , diluted to 1 
N H 2 S0 4 , and then hydrolysed at 121 °C (Fry, 2000). After neutraliza- 
tion with Ba(OH) 2 and removal of BaS0 4 , samples were subjected to 
methanolysis and trimethylsilylation followed by GC-flame ionization 
detection (FID) on a Hewlett Packard 5890. 

Bioinformatics 

MAPMAN and PAGEMAN software (Usadel et al, 2005, 2006) 
(http://mapman.gabipd.org) was used for statistical analysis and 
graphical representation of metabolic and signalling pathways. 
Dynamic analysis of the results described herein can be recapitu- 
lated by installing MAPMAN on a local computer and uploading 
the hlq data set in Supplementary Datafile 1 (sheet 1 ) at JXB online. 



Results 

The harlequin mutant manifests pleiotropic phenotypes 
affecting ABA- and auxin-inducible reporter gene 
expression, cell elongation, epidermal morphogenesis, 
and carbon partitioning into primary (starch, cell wall) 
and secondary metabolites (anthocyanins, lignin) 

The hlq mutant has an extreme dwarf phenotype, with pleio- 
tropic effects on root and shoot morphology, resulting in brittle 
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leaves and stems and more friable callus (Balasubramanian, 
2003), radial root tip swelling, anisotropic epidermal cell 
expansion, root hair defects, and abnormal cell collapse and 
associated callose accumulation (Supplementary Fig. SI at 
JXB online) (Subramanian et al, 2002). The shape of the cells 
in the root cortex and also the epidermis were irregular, unlike 
the regularly shaped spherical cells in the wild-type tissues 
(Fig. 1). The hlq hypocotyl and root epidermal cell lengths 
were drastically reduced and the columnar cell files were dis- 
rupted (Supplementary Table SI; Fig. S1H), suggesting that 
the dwarf phenotype of hlq is due to a reduction in cell size, 
and not in cell number, consistent with the results of others 
indicating that TopoVI is primarily involved in cell expansion 
by endoreduplication (Breuer et al, 2007; Sugimoto-Shirasu 
et al, 2005). The leaf pavement cells were less articulated 
than the degree of interdigitation seen in the parental type, 
with the presence of gaps between the pavement cells simi- 
lar to the gaps found in the epidermis of hlq hypocotyls and 
roots (Subramanian et al, 2002). Flower organs of hlq plants 
showed a rough epidermis compared with the smooth surface 
of the parental-type petals (Supplementary Fig. S1F). These 
phenotypes are consistent with a role for the HLQ gene prod- 
uct in a developmental^ modulated pathway that impinges 
on ABA, light, and cell wall signalling, similar to ABI8IKOBI 
ELD1, ProDc3:GUS expression was not specifically confined 
to either the root trichoblasts or atrichoblasts, as is seen for 
the WER or GL2 genes (Lee and Schiefelbein, 1999). GUS 
staining of median cross-sections of the root and hypocotyl of 
2-week-old hlq plants indicated that the ectopic proDc3:GUS 
expression also extended into the cortex (Fig. 1A). 

The hlq mutant is chlorotic, with only 50% of the chlo- 
rophyll content of the parental type when grown in light 




1 



Fig. 1. The hlq mutant has ectopic expression of the proDc3:GUS 
reporter gene. (A) Cross-section through the hypocotyl of 2-week-old hlq 
and parental-type (B) plants, stained for GUS. Parental-type plants do 
not stain for GUS after 1 0 d. hlq shows ectopic staining not only in the 
epidermal cells but also in the cortex. The arrows indicate collapsed cells 
in the epidermal cell file, and the cell shapes are also irregular. (C) Ectopic 
expression of proDc3:GUS in a 7-day-old hlq mutant primary root. (D) 
Wild-type expression of proDc3:GUS. Scale bar=500 |j.m. 



(Subramanian et al, 2002), suggesting that cellular metabo- 
lism may be broadly affected. When hlq pollen was germinated 
on artificial media, the pollen tubes grew at rates comparable 
with that of the wild type (data not shown). However, out- or 
back-crosses of hlq did not yield seeds, nor did manual pol- 
lination of hlq mutant flowers with the parental-type pollen. 
The hlq mutant plants failed to survive on soil or on basal 
minimal agar media, but under conditions of 0.8% sucrose 
supplementation, hlq homozygous plants could survive to the 
flowering stage, although vegetative tissue lost chlorophyll 
and leaves appeared brown/red in colour (data not shown). 

The ability of hlq mutants to survive to flowering when 
grown on sucrose prompted investigation of carbon parti- 
tioning into starch, since it serves as a short-term carbohy- 
drate reservoir. The hlq mutant in the Ler genetic background 
appeared to accumulate anthocyanins in the shoot meris- 
tematic region (Supplementary Fig. SIC at JXB online), 
reminiscent of the de-etiolatedl constitutive photomorpho- 
genesislfusca class of mutants and transgenic Arabidopsis 
overexpressing ABB, ABM, and ABI5 transcription fac- 
tors (Finkelstein et al, 2002), but not when in the abi2-l 
parental background (data not shown). This is possibly 
because the abi2 mutation has been shown to suppress sugar 
induction of the transcription factors PRODUCTION OF 
ANTHOCYANIN PIGMENT1(PAP1)IMYB75 and PAP2I 
MYB90 ~6- and 9-fold, respectively, relative to wild-type Ler 
(Luo et al, 2012). 

Sugar partitioning into starch was assessed by staining 
hlq mutants with IKI (Kurata and Yamamoto, 1998) after 
10 h of darkness. The hlq mutants showed patchy staining of 
starch granules in the hypocotyls and in the leaves (Fig. 2A), 
while the wild type was completely devoid of starch (Fig. 2B). 
The ectopic deposition and relative abundance of starch in 
shoots of hlq mutants was reminiscent of the proDc3:GUS 
staining pattern. However, counterstaining the GUS-stained 
hlq seedlings with IKI showed no obvious concordance 
between cells expressing GUS and those accumulating starch. 
Furthermore, the bulged and collapsed hypocotyls cells did 
not correlate with the presence of starch (Fig. 2C). Together, 
these results suggested that hlq has either decreased starch 
breakdown or increased starch synthesis rates, or defects in 
carbon metabolism or partitioning, which might also affect 
flux into secondary metabolites such as phenolics (e.g. antho- 
cyanins and lignins). Starch was subsequently quantified in 
hlq and two new alleles (see below). 

Cano-Delgado et al. (2000) showed that cell elongation 
defects lead to ectopic lignification. Their results were inde- 
pendently confirmed here with the korriganl mutant (data 
not shown) and ectopic lignin was tested for within the pro- 
scutel-HAtCES6 cellulose-deficient mutant (Fagard et al, 

2000) and boterol-1, a mutant allele of a katanin-like protein 
thought to be involved in microtubule stability (Bichet et al, 

2001) . It was found that both mutants displayed ectopic ligni- 
fication within the root (Supplementary Fig. S2D, E at JXB 
online). In order to investigate the relationship between the 
cell wall and hormone- and stress-regulated gene expression, 
the expression of the proDc3:GUS reporter was examined 
in the prcl-1 and botl-1 cell wall mutants. The proDc3:GUS 
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Fig. 2. Ectopic localization of starch in hlq mutants. Staining for starch 
granules by Kl shows ectopic localization of starch in the hypocotyl and 
leaves of hlq (A, white arrows), whereas there is no staining in the wild-type 
except at the meristem (B, white arrows). A close-up view of hlq hypocotyl 
(C) shows that the starch localization does not strictly correspond to hlq 
specific cell types— neither a bulged, radially expanded cell (large black 
arrowhead) nor a collapsed cell (small black arrowhead). 

marker gene was crossed into the prcl-1 and botl-1 mutant 
backgrounds using a Ler double insertion line (P series) (Chak 
et al, 2000) as one parent. In the F 2 progeny, the expres- 
sion patterns of proDc3:GUS in mutant individuals were 
visualized with the chromogenic substrate X-Gluc. Roots of 
prcl-1 had varying degrees of constitutive GUS expression 
(Supplementary Fig. S2B), but were clearly distinguishable 
from wild-type plants (Supplementary Fig. S2A) on the basis 
of their proDc3:GUS staining intensity. Individuals of botl- 
1 also displayed a range of proDc3:GUS expression pheno- 
types (Supplementary Fig. S2C). The extent of proDc3:GUS 
expression correlated with the extent of severity of the dwarf 
phenotypes in the botl-1 population (based on root length). 
These results taken together suggest that cell wall mutants 
phenocopy the ectopic expression of proDc3:GUS seen in hlq 
and raise the question of whether the hlq mutant also mani- 
fests ectopic lignification. 

Significantly, 2 weeks after germination, hlq mutant plants 
showed conspicuous swelling of the root elongation zone 
(Supplementary Fig. SID at JXB online), whereas swelling 
was less obvious earlier in development (data not shown). 
The swelling appeared similar to the phenotype of cellu- 
lose deficiency for radially swollenl, which encodes a glyco- 
syl transferase (Arioli et al, 1998), and the phenotype of 
the cytokinesis defectivel mutant, which encodes a GDP- 
mannose pyrophosphorylase required for A^-glycosylation 
(Lukowitz et al. , 2001). Treatment of wild-type seedlings with 
the herbicide dichlobenil (DCB; 2,6-dichlorobenzonitrile), 



an inhibitor of cellulose synthesis, phenocopied the swelling 
observed in the hlq root elongation zone. In the presence of 
1 uM DCB, wild-type plants developed a swollen root dif- 
ferentiation zone with few root hairs (Supplementary Fig. 
S3F). A similar effect was induced by treatment with 1 uM 
tunicamycin, which interferes with the synthesis of the core 
glycan chain attached to the Asn-Xaa-Ser/Thr motif of 
AMinked glycoproteins, which resulted in severe dwarfing, 
abnormal swelling of the root, and epidermal cell bulging 
(Supplementary Fig. S3A), as described by Lukowitz et al 
(2001). Root hairs emerged closer to the root tip, suggesting 
inhibition of cell elongation that led to a shorter root elon- 
gation zone. The hypocotyls of wild-type plants grown on 
tunicamycin were also short and radially swollen, with abnor- 
mally bulged epidermal cells (Supplementary Fig. S3B) simi- 
lar to those observed (Fig. 2C) for untreated hlq seedlings. 
Furthermore, when stained for proDc3:GUS activity, the 
tunicamycin-treated wild-type plants exhibited ectopic GUS 
expression (data not shown), consistent with results from cell 
wall mutants. 

Similar to the phenocopy of hlq for root swelling by 
inhibitors of cellulose biosynthesis, root tips of tunicamycin- 
treated seedlings had ectopic lignification, whereas root dif- 
ferentiated regions were normal (Supplementary Fig. S3C 
at JXB online; compare differentiated root region on the left 
with root tip on the right). Staining for lignin with phloro- 
glucinol showed ectopic staining in the cotyledons and roots 
of hlq (Fig. 3A, B), while the wild type showed only a faint 
staining in the vasculature (Fig. 3C). Histochemical staining 
of hlq mutants for the presence of suberin using Sudan Red 
did not show differences from the wild type. Finally, stain- 
ing with Ruthenium Red, which specifically identifies pectins 
(Shevell et al, 1994), indicated ectopic pectin accumulation 
throughout the hypocotyl and cotyledons of hlq mutants 
(Fig. 3E), which compares with the detection of pectin only 
in the vasculature, meristem, and leaf primordia of wild-type 
plants (Fig. 3D). Taken together, these results suggested that 
the hlq mutant has defects in carbon partitioning to starch 
and in cell wall biosynthesis impacting primary (cellulose, 
pectin) and secondary carbon metabolism, for example lignin 
and anthocyanins. 

hlq is an allele of BRASSINOSTEROID 
INSENSITIVE3/ROOT HAIRLESS3/HYPOCOTYL6/ 
TOPOISOMERASE6 SUBUNIT B 

The hlq mutant was previously mapped to the short arm of 
chromosome 3, ~6.6 cM south of the SSLP marker nga!62 
(Subramanian et al. , 2002). The locus was further fine-mapped 
by analysing another 3832 F 2 recombinant chromosomes 
segregating as homozygous mutants from two self-pollinated 
F[ plants (976 and 940 recombinant progeny, respectively) 
obtained from a cross between hlq/+ heterozygotes and the 
ecotype Col-0. The results placed hlq in an ~1 12 kbp interval 
within bacterial artificial chromosome (BAC) clones F3H11 
(>coordinate chr3:7213133), MOE17, and MDF22 ^coor- 
dinate chr3:7325482), a region that contains 30 annotated 
genes. The first indications as to the possible molecular defect 
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Fig. 3. Ectopic expression of lignin and pectin in hlq seedlings. Phloroglucinol staining for lignin in hlq shoot (A) and roots (B) shows ectopic deposition 
of patches of lignin (arrows). Wild-type (C) staining for lignin showed a faint signal only in the vasculature. Bars=1 mm for (A, C); 2 mm for (B). (D, E) 
Ruthenium Red staining for pectin. Ten-day-old wild-type seedlings (D) showed faint staining in the vasculature, meristem, and leaf primordium, whereas 
there was ectopic deposition of pectin in hlq (E) throughout the vasculature, leaves, and hypocotyl. Bar=1 mm. 



in hlq came from the observation that another mutant, hyp6, 
characterized by K. Sugimoto-Shirasu (personal commu- 
nication, XHIth International Conference on Arabidopsis 
Research) mapped near to hlq and shared similar phenotypes 
of anisotropic growth, radial swelling of epidermal cells, 
strongly reduced hypocotyl elongation, immature trichome 
development, suppressed root hair formation, and enriched 
cell wall pectins (Balasubramanian, 2003). 

A candidate G— >A mutation (TAIR10 coordinate chr3: 
7267232) that abolished the intron 7 donor consensus 
(nucleotide 1 of 300) in the pre-mRNA of At3g20780IBIN3l 
H YP6IRHL3I TOP6B was identified by Illumina HiSeq 2000 
whole-genome DNA resequencing of hlql+ heterozygotes. 
Supplementary Fig. S4 at JXB online shows the results of 
Sanger sequencing with a double peak at the mutation site in 
amplicons from hlql+ heterozygous template genomic DNA, 
confirming the Illumina short read analysis. The predicted 
effect of the mutation, assuming no splicing of intron 7 fol- 
lowed by translation of the aberrant 2512 nucleotide mRNA 
(based on TAIR10 annotation), is a truncated protein lacking 
472 C-terminal amino acid residues of the normal full-length 
670 amino acid protein, with the peptide IIIYSYQV added 
to Met 198 before a nonsense UGA stop codon at nucleotide 
25 of intron 7. The functional consequence of the translated 
mutant mRNA would be mutation of a highly conserved ali- 
phatic residue, two residues removed from the GXG motif 
of the ATP-binding site in the ATPase-like domain (amino 
acids 49-215; pfam accession PF02518; http://www.ncbi. 
nlm.nih.gov/Structure/cdd/cddsrv.cgi?uid=238030) and loss 
of the DNA TopoVI subunit 6B transducer domain (amino 
acids 398-557; pfam accession PF09239) (Corbett and 
Berger, 2006). 

In order to provide supporting evidence for the G— >A 
transition mutation being causal for the hlq phenotypes, the 
effect of the mutation was characterized at the mRNA level. 
Figure 4 shows the results of a total RNA gel blot from seed- 
lings segregating hlq/+ heterozygotes (2:1 ratio with wild- 
type plants) and homozygous hlq mutant plants, probed with 
a 1.5 kbp amplicon encompassing most of the At3g20780 
mRNA (2212 nucleotides). There was a reduction in the rela- 
tive abundance of the wild-type At3g20780 mRNA (-2300 



nucleotides calculated by linear regression of the mobilities 
of the RNA size markers) in heterozygotes, and a complete 
absence of the wild-type mRNA species in hlq homozygous 
plants. This was replaced by a diffuse band, ~300 nucleotides 
larger (Fig. 1A), which is consistent with the prediction that 
the mutation disrupts splicing of the 300 nucleotide intron 7 
in At3g20780 pre-mRNA and that the hlq mutation is a null 
allele. Figure 4B shows the result of an RT-PCR experiment 
using cDNA templates prepared from wild-type total RNA, 
hlq/+ heterozygote RNA, and hlq (-/-) homozygous mutant 
RNA. This further substantiates the previous results, as an 
amplicon of ~380bp corresponding to the 300 bp predicted 
intron 7 plus adjacent sequences was seen only in hlql+ het- 
erozygote and hlq homozygote (-/-) genotypes, and a smaller 
amplicon of ~ 80 bp was seen both in hlql+ heterozygotes and 
control templates plus reverse transcriptase (+RT). Because 
the size of the wild-type amplicon was larger than the pre- 
dicted 55 bp and only one cDNA clone (AJ297843) has been 
described (Hartung and Puchta, 2001), it is plausible that an 
alternative splice donor found 21 nucleotides downstream in 
intron 7 might be functional, which would add seven amino 
acid residues to the translated wild-type Top6B. 

Despite performing several backcrosses of the original 
mutant isolate to the parental line to remove extraneous 
mutations, because of close linkage, 12 other candidate het- 
erozygous SNPs were found in the mapped interval by whole- 
genome resequencing of hlql+ heterozygotes. This formally 
raises the possibility that the other mutations might contrib- 
ute to the pleiotropic hlq phenotypes, notwithstanding the 
fact that none of these SNPs falls within annotated exons or 
at intron donor-acceptor consensus sites. In order to provide 
conclusive evidence for a causal link between the hlq muta- 
tion and the pleiotropic phenotypes, two T-DNA insertion 
lines were obtained (SALK_140704 and SALK_024455) 
which disrupt exons 4 and 12 of At3g20780, respectively. 
Segregation of dwarf individuals with pleiotropic defects 
in root hair abundance and epidermal cell morphology 
was observed (Supplementary Figs S5, S6 at JXB online). 
Significantly, homozygous mutants of both these lines accu- 
mulated anthocyanins in meristems (Supplementary Fig. 
S5), as seen for hlq (Supplementary Fig. SIC). Based on the 
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hlq-l (G-»A, chr3.7267232) 

Fig. 4. Characterization of the At3g20780 mRNA in the hlq mutant demonstrates aberrant intron 7 splicing. (A) RNA gel blot of At3g20780 shows an 
aberrant mRNA in hlq homozygotes. A 20 |ig aliquot of total RNA from seedlings of wild types Col-0, Ler, and HLQ+/+ (afa/2-7 parental background), 
hlq/+ heterozygotes, and the homozygous hlq/hlq mutant were probed with a 1 .5 kbp fragment of the At3g20780 cDNA; the lower panel shows 
equal loadings by ethidium bromide staining of the denaturing agarose gel before blotting. Band sizes of the wild-type and hlq At3g20780 mRNAs 
are labelled on the left and right sides, respectively. (B) Agarose gel of RT-PCR amplicons from cDNA prepared from total RNAs of control (Col-0, Ler, 
and HLQ/HLQ), hlq/+ heterozygotes, and hlq/hlq homozygotes (-/-) shows retention of intron 7 (355 bp predicted) in hlq genotypes, and absence of 
intron 7 (-80 bp amplicon, 55 bp predicted) in controls. '+RT', with reverse transcriptase; '-RT reverse transcription step omitted from the amplification 
protocol. (C) Cartoon showing the exon-intron-untranslated region (UTR) structure of At3g20780, the position of primers (PF, PR; horizontal arrows) 
used for PCR in B, the position of the hlq-1 mutation (vertical arrow) at the intron 7 splice donor site, and positions of T-DNA insertions (triangles) in lines 
SALK_024455C (hlq-3) and SALK_1 40704 (hlq-2). 



nomenclature of Schiefelbein and Somerville (1990), hlq and 
the tested T-DNA insertion lines had few normal root hairs 
of at least 100 um in length, and sparse spike (non-branched) 
trichomes (Subramanian et al, 2002) were seen on primary 
leaves. The initiation of root hairs was normal from the basal 
end of the trichoblasts (towards the root apex), suggesting 
that the epidermal cells are not impaired in their ability to 
determine polarity. This observation is in contrast to the mid 
mutant of the TopoVI regulatory subunit, which has been 
reported not to express an apical polarity marker ROP2 asso- 
ciated with root hair initiation (Kirik et al, 2007). 

Qualitatively the root hairs of hlq (not shown) and the 
T-DNA insertion lines were observed to bulge or branch 
(Supplementary Fig. S6 at JXB online). The roots and shoots 
of the knockout mutants had a visually rough appearance, 
and the root differentiation zone and the distal region of the 
elongation zone were often swollen (Supplementary Figs S5D, 
S6). Homozygotes of both T-DNA insertion lines had very 
similar pleiotropic phenotypes to those of hlq mutants, but 
not as severe in terms of dwarfing of roots and shoots (Fig. 5). 
Interestingly, the presumed homozygous line SALK_024455C 
was found not to be homozygous, with only a few seeds out of 
~60 obtained from the stock centre showing the dwarf phe- 
notype. After verifying heterozygous lines, complementation 
tests were performed by crossing to hlql+ heterozygotes, and 




Fig. 5. T-DNA knockout lines SALK_1 40704 and SALK024455C identify 
two new hypomorphic alleles of hlq based on severity of shoot and root 
dwarfism and other pleiotropic effects on morphogenesis. Segregating 
mutants were plated on 0.8% sucrose minimal media and grown for 6 
weeks. Bar=1 cm. 



non-complementation (predicted 3:1 wild type:mutant segre- 
gation) was observed in F[ seedlings (Supplementary Fig. S7). 
PCR genotyping of these individual plants showed that they 
were heterozygous for the T-DNA insertions, which rules out 
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self-pollination of hlq/+ as an alternative interpretation of 
the results. % 2 analysis of the numbers of observed segregat- 
ing mutants in the Fj generation from crosses of heterozy- 
gous lines for the T-DNA insertions of SALK_024455C and 
SALK_140704 with hlql+ provided conclusive evidence that 
hlq is an allele of BIN3IH YL6IRHL3ITOP6B (Table 1). Since 
the T-DNA alleles have not been described previously, they 
are designated based on decreasing phenotypic dwarf severi- 
ties (Fig. 5) as hlq-2lsalk_140704 and hlq-3/salk_024455c, 
respectively (Fig. 4C). Analyses of absolute (Table 2) and 
relative (Fig. 6) starch accumulations in leaves and roots of 
the hlq-1, hlq-2, and hlq-3 mutants compared with isogenic 
controls validated prior staining results (Fig. 2A, C), and fur- 
ther showed that excess starch accumulation in mutants is not 
abated during the night, when starch was normally metabo- 
lized to give an ~7-fold drop in wild-type controls (Table 2). 
The strong correlations between fold starch accumulations in 
mutant leaves (Table 2) or roots (Fig. 6C) taken together with 
the observed severities of dwarf phenotypes (Fig. 5) supports 
the interpretation that hlq-1 is a null allele (Fig. 4A) and that 
hlq-2 and hlq-3 are somewhat leaky or hypomorphic alleles. 

Transcriptome profiling of the hlq mutant and meta- 
analysis of other TopoVI mutant data sets reveals 
pleiotropy at the levels of chromatin and signalling 
of light, ROS, hormone, sugar, and calcium affecting 
carbon metabolism 

Since the HL QIBIN3IH YL6IRHL3I TOP6B gene has pre- 
viously been described as regulating several molecular 



processes and signalling, it was endeavoured to take a systems 
approach to better understand the nature of HLQ effects, by 
analysing deregulated gene expression in hlq mutant seedlings 
using microarrays. The data are provided in Supplementary 
Datafile 1 at JXB online. Of the 25 673 genes on the micro- 
array, ~23% were significantly (/ > <0.05) misregulated in hlq 
compared with control, with 305 and 365 genes being down- 
regulated and up-regulated, respectively, >2-fold (average 
/ > <0.002). Using an arbitrary cut-off for fold effects of 1.74 
(log 2 =0.8), 514 genes were identified as down-regulated and 
559 up-regulated in hlq (average P<0.003). Supplementary 
Fig. S8 shows a 'volcano plot' for FC of all probes as a func- 
tion of statistical significance (P-values of FC). 

The Java desktop application 'Exploratory Gene 
Association Networks (EGAN) (Paquette and Tokuyasu, 
2010), which provides a contextual graph visualization of 
transcriptome results, was employed and an intriguing obser- 
vation was made: excluding 24 cases of gene duplications 
(where co-regulation is assumed to be due to duplicated regu- 
latory elements), there was a statistically significant chromo- 
somal adjacency for 27 sets/pairs of up-regulated and 24 sets/ 
pairs of down-regulated genes, or ~5% of the 1073 most dif- 
ferentially expressed genes (P<10"~ 21 , % 2 test using 1078 'con- 
trol' genes with average FC symmetrically distributed about 
zero; Supplementary Table S2 at JXB online). This number 
of gene pairs was more than three times the number of cases 
that would be predicted based on observed discordant pairs 
(14 pairs of adjacent differentially expressed genes that had 
opposite signs for FC, which fit well to a binomial model 
of chance occurrence). This observation suggested that the 



Table 1 . % 2 analysis of mutant phenotype segregation (non-complementation) observed in F 1 progeny of crosses between 
SALK_24455C/+ and SALKJ 40704/+ with hlq/+ heterozygous genotypes supports that HLQ Is an allele of BIN3/HYP6/RHL3ATOP6B 



Phenotype Observed SALKJ24455 3 Expected" Observed SALK_140704 C Expected" 

Mutant 9 12 47 48.75 

Wild type 39 36 148 146.25 

Total individuals 48 1 95 

X 2 P-value 0.32 0.77 



a Six separate crosses. 

b 1 mutant:3 wild type segregation null hypothesis for non-complementation, df=1 . 
c Eleven separate crosses. 

Table 2. Quantitation of starch (mg starch per g FW) in different organs and in response to light in a series of hlq mutant alleles 



Genotype 


Tissue, and time of harvest 






Ratio of leaf starch in 
mutant/wild type (fold 
effect) 


Leaf, end of day 


Leaf, end of night 


Root, end of day 


Root, end of night 


End of day End of night 


hlq- 1 (-/-) homozygote 


15.94 


13.56 


5.63 


4.41 


4.1 28.6 


HLQ (+/-) heterozygote/WT 


3.89 


0.47 


2.62 


1.74 




salk_140704/hlq-2 (-/-) homozygote 


17.18 


9.76 


2.65 


1.72 


2.8 12.6 


SALKJ 40704 (+/-) het/WT 


6.09 


0.78 


1.43 


1.10 




salk_25544c/hlq-3 (-/-) homozygote 


9.88 


9.81 


3.25 


2.28 


1.9 11.9 


SALK_25544C (+/-) het/WT 


5.14 


0.82 


1.58 


1.06 




Wild type average (+SEM) 


5.04 + 0.64 


0.69 + 0.11 


1.88 + 0.37 


1.30 + 0.22 




Mutant average(+SEM) 


14.33 + 2.25 


11.04 + 1.26 


3.84 + 0.91 


2.80 + 0.82 
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Fig. 6. Relative starch accumulation in leaves (A, B) and roots (C, D) of 
hlq mutant alleles compared with isogenic controls at the end of day (A, 
C) and end of night (B, D) diurnal cycles. Asterisks indicate significantly 
different from the wild type (*P<0.05; **P<0.01 , Student's paired f-test, 
equal variance assumed). Error bars are the SEM, n=4. 



TopoVI complex containing HLQ/TOP6B acts in part by 
chromatin remodelling of proximal clusters of two and three 
genes, as has been observed for the differentiation of haema- 
topoietic stem cells to erythroid and neutrophil cell types 
(Kosak et al. , 2007). Consistent with this interpretation is that 
56% of those immediately adjacent co-regulated genes are 
encoded on opposite strands, ruling out co-transcriptional 
control per se (e.g. by read-through of RNA polymerase II). 
The average distance between the distal ends of the co-reg- 
ulated genes was 8 kbp (Supplementary Table S2). Previous 



studies showed that the TopoIIa-like RHL1 and BIN4/MID 
proteins bind each other, TOP6A (but not TOP6B), and DNA 
(Breuer et al, 2007; Sugimoto-Shirasu et al, 2005). It is of 
note that Simkova et al. (2012) showed by chromatin immu- 
noprecipitation that Arabidopsis TopoVI subunits RHL1 and 
TOP6A directly bind to the promoter/transcription initiation 
sites of a pair of adjacent genes encoding unknown domains, 
Atlg24145 and Atlg24147, claiming that the TopoVI com- 
plex functions directly in initiation and elongation of tran- 
scription for different classes of genes subject to regulation by 
singlet and peroxide pathways of ROS. 

Given the recent claim that TOP6A/CAA39/AtSPOH-3/ 
RHL2/BIN5 is an integrator of singlet oxygen and H 2 0 2 
stress response pathways and therefore regulates expression 
of several thousand genes in Arabidopsis (Gadjev et al. , 2006), 
it was first sought to compare those reported '0 2 - and H 2 O r 
responsive genes subject to regulation by TOP6A/CAA39 
(Simkova et al, 2012) with first-generation transcriptome 
microarray studies on bin5ltop6a and bin3ltop6b (Yin et al, 
2002) as a framework for interpretation of the transcrip- 
tome results with hlqltop6b. Of the 255 genes claimed to be 
down-regulated in the caa39ltop6a mutant, 13 genes (or only 
5%) were strictly concordant with the 314 genes claimed to 
be down-regulated in bin5ltop6a (Supplementary Table S3 at 
JXB online) (but with the caveats that bin5 results were from 
a smaller number of genes interrogated by the Affymetrix 
8.3K GeneChip®, and caa39 is claimed to be a weak allele). 
This observation raised questions about interpretations of 
specificity drawn by previous microarray studies of TopoVI 
mutants. 

Transfer from dark to light of 5-day-old top6alcaa39 
mutants in the fluorescent in blue light (flu) mutant back- 
ground triggers photo-oxidative stress, resulting in mis- 
expression of 1093 genes (Simkova et al, 2012), similar to 
the number of significantly misregulated genes that were 
observed in hlq. This situation permits a reasonable analy- 
sis for agreement across experiments. Supplementary Table 
S3 at JXB online provides a comparison of genes displaying 
altered transcript levels for top6alcaa39 (Simkova et al. , 2012) 
with the present results for hlqltop6b. For the 5 1 genes (clus- 
ter 1) reported to be up-regulated in top6alcaa39 mutants 
unrelated to singlet oxygen signalling (and enriched for genes 
related to DNA repair), 13 were significantly up-regulated in 
hlq and one was significantly down-regulated; this is a very 
good agreement in terms of qualitative effect, but represents 
only about a quarter of the results that would be expected, 
assuming reproducibility. 

For the 113 genes claimed to be induced by '0 2 in a TopoVI- 
dependent manner (cluster 4), 22 genes were observed to be 
significantly up-regulated and 12 genes significantly down- 
regulated, for an overall concordance of 30% differential 
expression, but not strictly (the null hypothesis of repro- 
ducibility was down-regulation in hlq, but more genes were 
up-regulated). For the cluster of 55 genes hyperactivated by 
'0 2 in top6alcaa39 mutants (cluster 5, indicated as TopoVI- 
repressed), 12 up-regulated and two down-regulated genes 
were observed in hlq, a similar concordance (25%) with the 
TopoVI-specific effects described above. Similarly for cluster 
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2, in which were indicated 25 down-regulated genes in top6al 
caa39 not subject to singlet oxygen regulation, four down-reg- 
ulated and two up-regulated genes were observed. This gener- 
ally fits the null hypothesis, but with a low concordance (24%) 
as observed for the singlet oxygen-specific genes claimed to 
interact with top6alcaa39. For the 199 genes (clusters 6 and 
8) indicated as singlet oxygen inducible and up-regulated in 
top6alcaa39 (i.e. indicated to be TopoVI repressible), good 
concordance was observed with only five down-regulated 
genes and 90 up-regulated genes (48% of the total). 

Turning attention to further data sets, an analysis of ~5500 
genes in which 316 genes were reported as down-regulated in 
both bin3ltop6b and bin5ltop6a mutants (Yin et al, 2002) was 
examined. Supplementary Table S3 at JXB online provides the 
MqfWT log 2 FC results for genes identified in that list. Good 
concordance was observed, with 137 genes being significantly 
misregulated in the hlqltop6b mutant (43%), and very good 
reproducibility, with only 6% (20) of these genes being up-reg- 
ulated and 117 genes down-regulated in hlq, as hypothesized. 
Taken together, this meta-analysis supports the hypothesis 
that TopoVI mutants reproducibly alter expression of specific 
gene sets that can provide clues to the underlying cross-talk in 
gene networks controlling growth and development. 

Finally, those genes (clusters 3 and 7) reported as down- 
and up-regulated, respectively, by singlet oxygen produced 
after a shift to high light conditions in the flu background, 
but which were claimed not to interact with the top6alcaa39 
genotype (Simkova et al, 2012), were compared. In this case, 
the null hypothesis in the meta-analysis is non-concordance 
(low reproducibility) for the hlq/top6b-specific effects. Of the 
203 genes reported as down-regulated by a shift to high light, 
39 were significantly up-regulated in hlq, whereas 24 genes 
were significantly down-regulated (total concordance 31%), 
similar to or higher than the overall percentage seen for the 



observed concordance for TopoVI-specific effects (clusters 1, 
2,4,5,6, and 8). 

Similarly, for cluster 7 (365 genes up-regulated by light 
shift but independent of Wp6alcaa39), 132 significantly up- 
regulated and 24 down-regulated genes were observed in 
hlqltop6b, a very good agreement of 43% overall and strict 
reproducibility in terms of expected up-regulation. Taken 
together, these results showing high hlqltop6b concordance, 
vis-a-vis claimed top6alcaa39 non-specific effects, suggest 
that the TopoVI modes of action include pathways other 
than the oxidative stress response. Despite the relatively 
good agreement for differentially expressed genes observed 
between the present results and two independent experiments 
that reported similar numbers (several hundred) of TopoVI- 
dependent down-regulated genes, meta-analysis across all 
three microarray experiments found only seven consistently 
misregulated genes, and strict concordance was marginal 
(Table 3). This finding reinforces the notion that the pleio- 
tropic nature of TopoVI mutants defies easy descriptions of 
function, which raises questions of whether TopoVI is indeed 
a key integrator of the principal pathways of brassinoster- 
oids (Yin et al, 2002), cell cycle control (Hartung et al, 2002; 
Sugimoto-Shirasu et al, 2002), or singlet oxygen (Simkova 
et al, 2012) signalling as previously claimed. Since hlq came 
out of a screen for misregulated ABA-inducible gene expres- 
sion (Subramanian et al, 2002), TopoVI activity probably 
targets many different signalling pathways and processes, 
either directly or indirectly. 

Given the pleiotropic nature of hlq, MAPMAN software 
(Usadel et al, 2005) (http://mapman.gabipd.org) was used 
to analyse empirically the transcriptome data set in terms of 
metabolic and regulatory pathways affected in hlq mutants. 
Supplementary Fig. S9 and Supplementary Datafile 1 at JXB 
online provide, in the form of a relative heat map and lists, 



Table 3. Genes significantly misregulated in three independent microarray experiments* in TopoVI mutants 



AGI 


Annotation, gene ontology 
process/function 


top6alcaa39P 


hlq/top6b log 2 FC 
(this work) 


Concordance? 


At5g42650 


CYP74A/DELAYED DEHISCENCE2; allene 
oxide synthase 


Singlet induction TOP6A 
dependent 


-0.67 


Yes 


At5g47370 


HAT2, homeodomain transcription factor 
induced by auxin; root morphogenesis related 


Singlet repressed, TOP6A 
independent 


-0.29 


Yes 


At5g57560 c 


TOUCH4/cell wall modifying; rapidly induced 
by environmental stimuli 


TOP6A repressed, independent 
of singlet induction 


1.29 


Not for bin3 


At1g01120 


KCS1 , 3-ketoacyl-CoA synthase, critical for 
fatty acid elongation in wax biosynthesis 


Singlet induced, TOP6A 
independent 


-0.49 


Not for singlet 


At2g44940 


DREB subfamily A-4 of Ethylene Response 
Factor/APETELA2 domain transcription factor 


Singlet induced, TOP6A 
independent 


-0.49 


Not for singlet 


At4g20780 


CALMODULIN LIKE 42, EF hand domain 
calcium and protein binding; involved in 
trichome branching, abiotic stress responses 


Singlet induced, TOP6A 
independent 


-0.27 


Not for singlet 


At2g1 7880 


Chaperone DnaJ heat shock protein, 
response to sugar 


Singlet repressed, TOP6A 
independent 


0.25 


No 



a Originally reported (Yin et al., 2002) as down-regulated >2-fold in both bin3/top6b and bin5/top6a mutants. 

b Genes classified (Simkova etal., 2012) as 'topo6a/caa39 induced' (cluster 1); 'topo6a repressed' (cluster 2); 'singlet repressed, TOP6A 
independent' (cluster 3); 'singlet induced TOP6A dependent' (cluster 4); 'singlet induced TOP6A repressed' (cluster 5); 'TOP6A repressed, 
independent of singlet induction' (clusters 6, 8); 'singlet induced, TOP6A independent' (cluster 7). 
c Reported (Yin etal., 2002) as down-regulated; brassinosteroid induction dependent on BIN3ATOP6B and BINSATOPOA. 
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respectively, the number and degree of 2575 metabolism- 
related genes that are significantly misregulated in hlq (aver- 
age P<0.005 for visibly coloured genes in Supplementary Fig. 
S9). Genes for photosynthesis, photorespiration, terpenoid 
biosynthesis, nitrogen assimilation, tetrapyrrole/chlorophyll/ 
haem/phytochrome chromophore biosynthesis, the reductive 
pentose phosphate Calvin/Benson/Bassham cycle, and starch 
biosynthesis were uniformly down-regulated in hlq seedlings. 
This generally correlated with observed chlorosis in the mutant 
(Fig. 5), but was inversely correlated with starch accumulation 
(Figs 2, 6), an unexpected finding that suggests feedback regu- 
lation of carbon metabolism. Supplementary Datafile 2 docu- 
ments the findings of specific genes, pathways, and processes 
significantly altered at the transcriptome level in hlq. 

It was of particular interest that one of the down-regulated 
(2.38 FC, PO.OOOOl) tetrapyrrole genes is At5gl3630l 
GENOMES UNCO UPLED5IMg- CHELA TA SE-H, encoding 



an enzyme of chlorophyll biosynthesis variously described as 
involved in plastid to nucleus retrograde signalling (Mochizuki 
et al, 2001; Susek et al, 1993) and a mediator of ABA sig- 
nal transduction, possibly as an ABA receptor modulating 
expression of WRKY-domain transcriptional repressors, 
ROS homeostasis, and lipid (3-oxidation (Shen et al, 2006; 
Wu et al, 2009; Shang et al, 2010; Jiang et ah, 2011; Tsuzuki 
et al, 2011; Xu et al, 2012). Remarkably, the downstream 
targets of GUN5, At2g33150l3-ketoacyl-CoA thiolasel KAT2 
and Atlg80840IWRKY40, are implicated in ABA signal- 
ling (Shang et al, 2010; Jiang et al, 2011; Liu et al, 2012), 
and both are significantly up-regulated in hlq mutants (2.96 
FC, /><0.0001; and 1.96 FC, P<0.002, respectively (Table 4). 
KAT2 is categorized in Supplementary Fig. S9 at JXB online 
as involved in branched-chain amino acid catabolism, which 
was over-represented with significantly down-regulated bio- 
synthetic genes and up-regulated catabolic genes (_P=0.003: 



Table 4. ABA biosynthesis and signalling genes with significantly altered expression in hlq mutant seedlings 



AGI 


Gene name 


Annotation 


Wg/WT log 2 FC 


P-value 


At1g15520 


ABCG40 


ABC transporter; ABA importer 


1.91 


0.0001 


At2g47130 


SDP.3 


Short-chain dehydrogenase/reductase ABA2-like 


1.73 


0.0001 


At2g33150 


KAT2/PED1 


Peroxisomal 3-ketoacyl-CoA thiolase, fatty acid 


1.56 


0.0001 






p-oxidation; positive ABA effector 






At1 g80840 


WRKY40 


Pathogen-induced transcription factor; ABA repressor 


0.97 


0.001 


At5g67300 


MYB44 


R2R3 MYB transcription factor; ABA sensitivity 


0.85 


0.01 


At1g56070 


LOS1 


Low response to Osmotic Stressl ; translation 


0.52 


0.001 






elongation factor 2-like; cold stress response 






At5g58670 


PLC1 


Phospholipase C1 ; ABA, drought, salt, cold response 


0.45 


0.05 


At4g26080 


ABI1 


ABA insensitivel ; protein phosphatase 2C 


0.35 


0.008 


At3g57530 


CPK32 


Calcium-dependent protein kinase32; phosphorylates 


0.31 


0.01 






ABA RESPONSE FACTOR4 






At1g33560 


ADP.1 


ACTIVATED DISEASE RESISTANCE1 ; NBS-LRR; 


0.31 


0.03 






interacts with ABI1 in drought response 






At1 g69260 


AFP1 


ABI5-Binding Protein; domain unknown function 


0.28 


0.03 


At3g24650 


ABI3 


ABA insensitive3; B3 domain transcription factor 


0.26 


0.02 


At2g40220 


ABI4 


ABA insensitive4; AP2 domain transcription factor 


0.23 


0.03 


At4g01026 


PYL7 


Pyrabactin-I_ike7; ABA receptor 


0.22 


0.04 


At5g45870 


PYL12 


Pyrabactin-I_ike12; ABA receptor 


0.22 


0.05 


At3g2501 0 


P.LP41 


Receptor Like Protein41 ; ABA hypersensitive to 


-0.21 


0.04 






chlorosis 






At2g32860 


LOS 15 


AtBGL2; Beta glucosidase33; ABA-GE hydrolase 


-0.24 


0.02 


At1 g80080 


TMM/RPL 1 7 


Too Many Mouths; stomatal development; 


-0.26 


0.02 






ABA insensitive to chlorosis 






At1 g3551 5 


HOS10/MYB8 


High response to Osmotic Stressl 0; ABA 


-0.33 


0.02 






hypersensitive 






At5g45340 


CYP707A3 


ABA 8'-hydroxylase; phaseic acid synthesis 


-0.46 


0.03 


At1g73000 


PYL3 


Pyrabactin-Like3; ABA receptor 


-0.46 


0.02 


At4g1 7870 


PYP.1 


Pyrabactin Resistentl ; ABA receptor 


-0.58 


0.01 


At5g53160 


PYL8 


Pyrabactin-Like8; ABA receptor 


-0.63 


0.0008 


At3g63520 


NCED1 


Nine-c/s-epoxycarotenoid dioxygenase 


-0.75 


0.0005 


At5g67030 


ABA1 


Zeaxanthin epoxidase; mutant ABA deficient 


-0.84 


0.002 


At4g15560 


CM 7 


1 -Deoxyxylulose 5-phosphate synthase; 


-0.87 


0.0002 






chloroplastos alteradosl ; ABA deficient 






At1g52400 


AtBG1 


Beta glucosidase18; ABA-GE hydrolase 


-1.01 


0.0002 


At2g41070 


DPBF4/EEL 


Dc3-Pro Binding Fctr4; bZIP; Enhanced Em Level 


-1.24 


0.0002 


At1 g64670 


CED1 


9-CIS EPOXYCAROTENOID DIOXYGENASE 


-1.60 


0.0008 






DEFECTIVE; BODYGUARD1 ; epidermal wax biosyn- 










thesis; alpha-beta hydrolase 
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Supplementary Datafile 1). Interestingly, several of the 
WRKY genes up-regulated in hlq (Supplementary Table S4) 
are also expressed specifically in trichomes ( WRKY8, -15, -18, 
and -33) (Jakoby et al, 2008) and involved in hormone sig- 
nalling cross-talk and oxidative stress response (WRKY25) 
(Hu and Ma, 2006). Similar to mutants that do not produce 
mature trichomes (Marks et al, 2009), many abiotic stress 
effectors were misregulated in hlq (Supplementary Datafile 1), 
including genes for ABA biosynthesis and metabolism, mem- 
bers of the PYRABACTIN-RESISTANT1 1 REGULATORY 
COMPONENTS OF ABA RECEPTOR11IRCAR11 class of 
ABA receptors, and genes with functional evidence as effec- 
tors of ABA signalling and biosynthesis (Table 4). 

Supporting the observed phenotypes of lignin and antho- 
cyanin accumulation in hlq (Fig. 3A, B; Supplementary Figs 
SIC, S5D at JXB online), many genes involved in aromatic 
amino acid and phenylpropanoid/lignin biosynthesis were 
up-regulated in hlq, whereas flavonoid biosynthetic genes 
were both up- and down-regulated (Supplementary Fig. S9; 
Supplementary Datafile 1), some of which are mono-oxyge- 
nases that have haem cofactors (see above). The strong up- 
regulation of NITRILASE2, which functions in an alternative 
pathway of wound- and pathogen-induced auxin (IAA) bio- 
synthesis from indole-3-acetonitrile (Huh et al, 2012) and 
degradation of tryptophan-derived glucosinolates, indicates 
that auxin or tryptophan amino acid homeostasis may be 
affected in hlq, or that the hlq defect in trichome development 
may be responsible for the observed (Supplementary Fig. S9) 



misregulation of glucosinolate biosynthetic genes (Jakoby 
et al, 2008). Since trichome-specific expression of the cell 
cycle checkpoint CYCLINB1;2 results in rescue of the top6al 
rhl2 and mid mutant trichome defects but not of their dwarf 
phenotypes (Kink et al. , 2007), it is likely that the misregula- 
tion of glucosinolate pathways in hlq mutants is a direct effect 
of TopoVI on carbon partitioning processes. 

With regards to the pleiotropic cell wall phenotypes of hlq, several 
genesencodingproteins(EXPANSrNS,ARABrNOGALACTAN 
PROTEINS) that modify or contribute to wall structure and/or 
metabolism were up-regulated (Table 5), providing correlative evi- 
dence for the functions of those genes in the cell wall phenotypes 
of hlq, including the increased abundance of pectins (Fig. 3E) and 
radial swelling of epidermal cells (Fig. 2C), which was phenocop- 
ied by the iV-glycosylation inhibitor tunicamycin (Supplementary 
Fig. S3 A, B at JXB online). Arabinogalactan proteins are involved 
in root and stem expansion (Willats and Knox, 1996; Ding and 
Zhu, 1997; Park et al. , 2003) and broadly implicated in plant devel- 
opment including programmed cell death (Gao and Showalter, 
1999) similar to the observed hlq phenotype, and possibly in ABA- 
inducible gene expression (Desikan et al, 1999). 

Characterization of hlq whole cell wall composition 
by infrared spectroscopy and quantitation of glycosyl 
composition 

The metabolism of callose ((3-1,3-glucan) within plant cell 
walls is essential to many developmental, physiological, and 



Table 5. List of cell wall genes most up- and down-regulated in hlq seedlings by transcriptome profiling 



AG I 



Annotation 



hlq/WT log 2 FC 



P-value 



Up-regulated wall biosynthesis 

At2g 18660 

At3g45970 

At2g22470 

At1 g35230 

At4g09030 

At5g6431 0 

At2g45220 

At3g48580 

At5g57560 

At4g23990 

At3g02230 

At4g2581 0 

At3g07160 

Up-regulated wall eatabolism 

At3g57510 

At4g13210 

At4g30270 

At5g20950 

Down-regulated wall metabolism 

At2g06850 

At4g37800 

At5g39310 

At2g37640 

At4g28250 



EXPANSIN-LIKE B3 
EXPANSIN-LIKE A2 
ARABINOGALACTAN PROTEIN 2 
ARABINOGALACTAN PROTEIN 5 
ARABINOGALACTAN PROTEIN 10 
ARABINOGALACTAN PROTEIN 1 
PECTINESTERASE 
ENDO-XYLOGUCAN TRANSFERASE 
XYLOGUCAN TRANSFERASE/TOUCH4 
CELLULOSE SYNTHASE/CSLG3 
UDP-I-Ara MUTASE/RGP1 a 

XYLOGUCAN ENDO-TRANSGLYCOSYLASE6/XTR6 
GSL10; GLUCAN SYNTHASE-LIKE1 0, p-1 ,3 CALLOSE 

POLYGALACTURONASE/ADPG1 
PECTATE LYASE 

MANNAN-XYLOSE HYDROLASE/MERISTEM5B 
(5-1 ,4-GLUCANASE 

ENDO-XYLOGUCAN TRANSFERASE/EXGT-A1 
ENDO-XYLOGUCAN TRANSFERASE 
EXPANSIN 24 
EXPANSIN A3 
EXPANSIN B3 



3.32 
1.29 
2.21 
1.14 
1.07 
1.00 
2.73 
1.20 
1.29 
0.84 
0.94 
1.06 
0.35 

1.42 
1.06 
1.07 
0.90 

-1.43 
-1.20 
-0.92 
-0.73 
-0.72 



<0.0002 
<0.0007 
<0.0002 
<0.0007 
<0.0003 

0.002 
<0.0009 
<0.0003 
<0.0002 
<0.0006 
<0.0006 
<0.0004 
<0.006 

<0.02 
<0.0003 
<0.0003 
<0.0008 

<0.0006 

<0.0004 

<0.002 

<0.002 

<0.003 



a Annotated in Rautengarten ef al. (201 1). RGP1 , REVERSIBLY GLYCOSYLATED POLYPEPTIDE1 
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stress-related processes (Chen and Kim, 2009). Production 
of callose is induced by mechanical stress and wounding, 
and under conditions of cellulose deficiency (e.g. in cell wall 
mutants). Some dwarf mutants ectopically deposit callose 
(Lukowitz et al, 2001), pectin (His et al, 2001), and the phe- 
nolic polymers suberin (Cheng et al, 2000) and lignin (Burk 
et al , 200 1 ), thereby compromising the stability and flexibility 
of the cell wall. Treatment of wild-type plants with tunicamy- 
cin resulted in the appearance of patches of callose accumu- 
lation within the root elongation zone (Supplementary Fig. 
S3D at JXB online), similar to the phenotype of hlq mutants 
(Subramanian et al, 2002). A comparison of callose content 
in hlq and the wild type during growth on sucrose-supple- 
mented agar media indicated that the hlq mutants showed 
a steady increase in callose accumulation over 15 d, with 
mutants accumulating up to five times the callose content of 
wild-type plants (Fig. 7). 

Yariv reagent binds arabinogalactan proteins and triggers 
accumulation of callose by induction of the callose synthase 
gene At3g07160/GLUCAN-SYNTHASE-LIKE10 and the 
transcription factor gene WRKY40 (Guan and Nothnagel, 
2004), which were both significantly up-regulated in hlq 
(Tables 4, 5). Similarly, the transcription factor gene MYB44 
was up-regulated in hlq (Table 4) and has been shown to regu- 
late callose production in Arabidopsis (Lu et al, 2011). Since 
the present results suggested that hlq mutants have a gener- 
alized perturbation in cell wall composition, hlq mutant cell 
walls were characterized by spectroscopic and biochemical 
methods. 

Along with soluble expansins and structural proteins 
including hydroxyproline-rich extensins and arabinoga- 
lactan proteins, four primary and secondary polymer classes 
make up the plant cell wall: cellulose, hemicellulose, pectin, 
and lignin. FTIR microspectroscopy, which measures the 
energy of asymmetric molecular bond vibrations, is a rapid, 
non-invasive method for detecting a range of polysaccha- 
rides affecting cell wall architecture in muro, based on their 




Fig. 7. Callose content of wild-type versus hlq during growth of seedlings. 
Mutants of hlq had five times more callose than the wild type after 1 5 
d of growth. Error bars represent ±SEM of three replicates. The results 
presented are for the abi2-1/hlq double mutant; however, the phenotypes 
described were independently verified on Ler/hlq single mutants. 



functional groups, including, but not limited to, carboxylic 
esters, phenolic esters, protein amides, and carboxylic acids. 
A digital subtraction plot of the hlq FTIR spectrum minus 
that of the wild-type spectrum from extracted cell walls is 
provided in Fig. 8A. The frequency-structure peak at wave- 
number 1727 cm 1 corresponds to the C=0 stretch of esteri- 
fied pectins, consistent with the qualitative result of pectin 
staining showing elevated pectins in hlq leaves (Fig. 3E) 
(Alonso-Simon et al, 2011). Secondary amide bonds attrib- 
uted to cell wall peptides (McCann et al, 1992) were observed 
at 1643 cm 1 corresponding to the C=0 stretch, at 1530cm 1 
for the N-H amide, and at 1229 cm" 1 for the C-N stretch. The 
peaks observed at 1138cm 1 and 1084cm 1 are probably due 
to the C-O-C glycosidic bonds of various polysaccharides 
that differ between hlq and the wild type. Importantly, the 
negative peak at 1055 cm 1 probably corresponds to the ring 
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Fig. 8. Fourier transform infrared (FTIR) microspectroscopy analysis of 
dehydrated cell wall extracts of hlq and the wild type. (A) Digital subtraction 
frequency-structure correlation chart. Difference (A) spectrum of hlq minus 
wild-type FTIR spectra. The peak at wavenumber 1727 cm" 1 corresponds 
to the C=0 stretch of esterified pectins; peaks at 1643cm" 1 (C=0 
stretch), 1530cm" 1 (N-H amide), and 1229cm" 1 (C-N stretch) correspond 
to secondary amide bonds attributed to cell wall peptides, or lignin for 
1643cm" 1 (McCann ef a/., 1992); peaks at 1 138cm" 1 and 1084cm" 1 are 
probably due to the C-O-C glycosidic bonds of various polysaccharides 
which differ between hlq and the wild type; the negative peak at 1055 cm" 1 
is the ring C-0 of cellulose, and 1010 cm" 1 is probably the ring C-C. (B) 
Principal component analysis of FTIR spectra separates the hlq and 
controls (WT and heterozygote samples) into two non-overlapping clusters 
which account for 98% of the variation. Inset: plot of the PC1 loadings 
versus wavenumber, showing that wavenumbers 1735, 1657, 1546, and 
1244cm" 1 (esterified pectins, lignin/peptides) are dominants in addition to 
cellulose at 1063cm" 1 (arrows). 



4232 I Mittalefa/. 



C-O of cellulose (Alonso-Simon et al, 2011) and homoga- 
lacturonan, whereas the positive peak at 1010 cm 1 probably 
represents the ring C-C stretch of various polysaccharides 
(pectins). These results suggest that there is a deficiency of 
cellulose or hemicelluloses and an increased abundance of 
pectins in the hlq mutant as compared with the wild type. 

The results of PCA discriminatory analysis on the hlq and 
wild-type wall spectra are shown in Fig. 8B. Most of the 
variation (98%) in FTIR spectral differences between cell 
wall extracts of hlq, hlql+ heterozygotes, and wild-type walls 
could be attributed to two PCs, with replicate hlq homozy- 
gous samples clustering separately in the lower left quadrant 
from hlq/+ heterozygotes, and wild-type samples cluster- 
ing broadly together in the other three quadrants (Fig. 8B). 
Plotting the PC that accounted for 95% of the variation as a 
function of wavenumber (Fig. 8B, inset arrows) showed peaks 
at 1735cm 1 (esterified pectins), 1657cm 1 (lignin), 1546cm 1 
(secondary peptide amide), 1244cm" 1 (C-N peptide amide), 
and 1063 cm 1 (cellulose or homogalacturonan), correspond- 
ing well to previous reports (McCann et al., 1992). It is con- 
cluded from this analysis that the hlq mutant is deficient in 
cellulose or hemicelluloses, or both, and has elevated levels 
of pectins, lignins, and cell wall proteins. This conclusion 
was further supported by the microarray results for cell wall 
structural proteins EXPANSIN-LIKE and arabinogalactan 
proteins (Table 5). 

In order to differentiate between the effects of hlq on cel- 
lulose and hemicellulose levels, the glycosyl composition of 
whole cell wall extracts subjected to methanolysis and tri- 
methylsilylation combined with GC with or without prior 




Wall Hydrolyzed Sugars 

Fig. 9. Glycosyl composition (mol%) of the whole cell wall fraction from 
wild type (WT) and hlq mutant Arabidopsis seedlings, without and with 
prior H 2 S0 4 swelling and hydrolysis (SAH) to differentiate cellulose and 
tightly bound mannosyl polymers from extractable (e.g. hemicellulose) 
polysaccharides. Cell walls of hlq had significantly decreased cellulose and 
compensatory increases in arabinosyl, rhamnosyl, and galactosyl residues 
typical of pectic arabinogalactans. Sugar residues were measured as TMS 
methylglycosides by GC-FID after methanolysis of whole walls (before 
and after SAH). Glc-SAH and Man-SAH represent the additional residues 
detected due to SAH. Error bars indicate ±SEM (n=2). Asterisks (*) indicate 
significantly different from both wild-type and hlq/+ heterozygote samples, 
P<0.02 (Student's two-sided r-test, equal variance assumed). H indicates 
significantly different from wild-type samples, P<0.02. 



swelling and hydrolysis of cellulose (SAH) using concentrated 
sulphuric acid was analysed. On a mole percentage basis of 
constituent sugar residues, the whole cell wall fraction of 
hlq seedlings contained significantly more abundant compo- 
nents of pectic arabinogalactans than the wild type, in par- 
ticular arabinose (5.5 mol% increase), rhamnose (1.3 mol% 
increase), and galacturonic acid (4.6 mol% increase), whereas 
hlq mutant walls contained significantly less cellulosic glucose 
(-14.5 mol% decrease) than those of the wild type (Fig. 9). 
This supports the finding that hlq accumulates more pectin 
in leaves than the wild type (Fig. 3E) and has a wall archi- 
tecture deficient in cellulose (Fig. 8), and the conclusion that 
cellulose deficiency and/or pectin overabundance resulted in 
severely compromised epidermal cell wall structural integrity 
(Fig. 2C; Supplementary Figs S1H, S6 at JXB online). 



Discussion 

Although previous claims (Hartung et al, 2002; Sugimoto- 
Shirasu et al, 2002; Schrader et al, 2013) on the role of 
TopoVI, found only in plants and archaebacteria, in the 
decatenation of endoreplicated chromosomes during cell 
growth are not under debate, based on the results obtained 
on hlq mutants, a more complex role for TopoVI in transduc- 
ing at the level of chromatin environmental and internal sig- 
nals controlling growth and development is proposed. Others 
have also recognized this possibility (Sugimoto-Shirasu 
et al, 2002; Yin et al, 2002; Breuer et al, 2007; Kirik et al, 
2007). The pleiotropic phenotypes of hlq were previously 
characterized, in particular that mutants are skotomorpho- 
genic, with deformed epidermal cells which stain for callose, 
abnormal and reduced root hairs and leaf trichomes, and sto- 
chastic ectopic expression of an ABA- and auxin-inducible 
transgene. Here hlq mutants were further characterized at 
the physiological and transcriptome levels as accumulating 
starch, callose, and anthocyanins, with cell walls abundant 
in pectins, lignin, and cell wall proteins. Moreover, despite 
the absence of endoreplication in Arabidopsis inflorescences 
(Galbraith et al, 1991), it was shown here (Supplementary 
Fig. SIB, F at JXB online), as has been shown by others 
(Sugimoto-Shirasu et al, 2002), that hlqltop6b mutants can 
produce mostly normal flowers but neither homozygous male 
nor female hlq sporophyte organs functioned properly in 
reproduction, notwithstanding artificial nutrient supplemen- 
tation or the previous claims that the hyp6ltop6b, rhl2ltop6al 
atspoll-3, and bin3-l alleles are null yet produce viable seeds 
(Sugimoto-Shirasu et al, 2002; Yin et al, 2002). Consistent 
with the notion that TopoVI is important for reproduction is 
the observation that in rice the TopoVI subunits are expressed 
maximally in pre-pollinated flowers (Jain et al, 2006). Here 
hlq and two additional alleles that vary in their phenotypic 
severity were isolated and it was shown that they are alleles of 
bin3/hyp6/rhl3 /Topoisomera&e 6B. 

The rhl3 allele of TOP6B is documented in TAIR (http:// 
www.arabidopsis.org) as a mutation of the splice acceptor in 
intron 17, which is stated to result in the transcription of an 
mRNA with at least three altered sizes that create a premature 
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stop codon in the last exon. Hartung et al. (2002) described 
an unnamed T-DNA mutant allele that disrupted exon 12 
and deleted 268 bp of the TOP6B gene including exonl3. On 
the basis of root calli death 3 weeks after growth induction, 
they argued that deficiency of AtTOP6B results in a general 
growth defect. Because of the absence, in the present case, of 
a correctly spliced intron 7 mRNA (Fig. 4) and based on the 
prediction that the hlq mutation sits adjacent to a highly con- 
served aliphatic residue two residues distal to the GXG motif 
of the ATP binding site in the ATPase-like domain (Corbett 
and Berger, 2006), it is argued here the hlq-1 mutation is a 
bona fide null allele. Phenotypic analysis of two new hypo- 
morophic hlq alleles (Fig. 5) supports this contention. The 
various descriptions by others of presumed null alleles dem- 
onstrate that, due to the pleiotropic nature of top6B mutants, 
limited phenotypic characterizations can give misleading 
interpretations of gene function. This notion is supported by 
the contradictory reports of bin3lhyp6ltop6b mutant effects 
on chromosome ploidy levels for several cell types (Hartung 
et al, 2002; Sugimoto-Shirasu et al, 2002; Yin et al, 2002) 
and the broad lack of concordance from meta-analysis of the 
hlq transcriptome with other previously published TopoVI 
mutant transcriptomes (Table 3; Supplementary Table S3 at 
JXB online), including other claimed null alleles. 

The Arabidopsis genome encodes a second type II topoi- 
somerase (At3g23890) that is differentially regulated from 
TopoVI (Breuer et al, 2007) and probably functions in mito- 
sis, T-DNA integration, and cell cycle regulation (Xie and 
Lam, 1994; Makarevitch and Somers, 2006), and may have 
some overlapping functions with TopoVI (Sugimoto-Shirasu 
et al, 2005). The claims that BIN4/MID, a DNA-binding 
component of TopoVI, functions primarily as a DNA dam- 
age response effector during post-mitotic endocycles (Breuer 
et al. , 2007) and transcriptional silencing (Kirik et al. , 2007) 
mediated through the cell cycle checkpoint effector Ataxia 
Telangiectasia-mutated and Rad3-related (At5g40820/ATR) 
and ectopically expressed markers At2g31320lpoly (ADP- 
ribose) polymer asel IP ARB 1 , At4g02390IPARP2, A T4G2 10701 
BRCA1, At5g20850IRAD51, and AT4G37490ICYCLINB1;1 
was not supported by the present transcriptome analysis of 
hlq, as none of these genes was altered (see Supplementary 
Datafile 1). Conclusions drawn from a meta-analysis of dif- 
ferent microarray data sets are not absolute. It cannot be 
ruled out that genotype-, tissue-, or methodology-specific dif- 
ferences between the present experiments and those of others 
contributed to lack of agreement at the transcriptome level. 
However, it is also noted that contradictory results have been 
reported for the role of BIN4/MID in transcriptional silenc- 
ing of heterochromatin (Breuer et al. , 2007; Kirik et al. , 2007). 
Since BIN4/MID does not physically interact with TOP6B/ 
HLQ (Breuer et al, 2007), the lack of similar molecular phe- 
notypes for DNA damage markers in hlq mutants suggests 
that alternative mechanisms exist for regulation of TopoVI 
activities independent of BIN4/MID. 

In this context, it can be recalled that hlq mutants were 
isolated on the basis of ectopic expression of the ABA- and 
auxin-inducible proDc3:GUS reporter gene in the epidermal 
cells of roots and hypocotyls (Subramanian et al, 2002). 



Previous results showed that only a subset of hlq root epi- 
dermal cells stained strongly with the viable stain fluores- 
cein diacetate (FDA). Remarkably proDc3:GUS staining 
exactly coincided with the FDA staining patterns, suggest- 
ing that FDA-marked epidermal cells specifically express 
proDc3:GUS and that these cells are metabolically hyper- 
active. The brassinolide insensitivity, skotomorphogenesis/ 
hypocotyl elongation defect, and root hairless phenotypes of 
the hlq mutant were previously described (Subramanian et al. , 
2002), consistent with the descriptions of other independent 
mutant allele phenotypes. The fundamental question is thus 
framed as 'chicken versus egg': is the hlqlbin3lhyp6lrhl3ltop6B 
mutation causal for the expression of proDc3:GUS, which 
would thus serve as a bona fide marker for Top6B function in 
gene regulation, or rather is proDc3:GUS expression a down- 
stream consequence of indirect internal stresses imposed by 
cell viability, wall defects/callose synthesis, or other second- 
ary processes disrupted in the hlq mutant? The present results 
showing altered biosynthesis of primary metabolites such as 
starch (Figs 2, 6; Table 2) and cell walls (Figs 3, 7-9) in hlq, 
as well as proDc3:GUS expression in cell wall mutants of 
different classes (Supplementary Fig. S2 at JXB online) and 
following inhibition of peptidoglycan and cellulose biosyn- 
thesis (Supplementary Fig. S3) are consistent with a primary 
(direct) effect of hlq on gene expression. It is also possible 
that the lignin accumulation phenotype (Fig. 3) could cause 
alterations in cell-cell communication and thereby disturb 
hormonal signals or sensitivities. 

The present meta-analysis of other transcriptome results 
showing the concordance of bin5ltop6a and caa39ltop6a with 
the hlqltop6b results further supports the chromatin remod- 
elling model, because numerous differentially expressed 
unrelated adjacent genes could be identified (Supplementary 
Datafile 1 at JXB online) such as At4g37250ILRR-NBS 
receptor kinase and At4g37260IMYB73; At5g23340IF-box 
and At5g23350/GRAM domain ABA responsive; Atlg24140l 
matrixin metallopeptidase and Atlg24150lformin homologue4; 
Atlg28370IERF-AP2 domainll and Atlg28380lnecrotic 
spotted lesionsl; ^4?/g<560S0/unknown and Atlg66090ILRR- 
NBS receptor kinase; Atlg74930IDREB subfamily A-5 and 
Atlg74940/DUF581; At2gl8690/unknown and At2gl8700l 
TSP 11;AT3G2 56001 EF-hand family protein and AT3G25610I 
haloacid dehalogenase; AT4G27652 and AT4G27654 (unre- 
lated unknowns); At4g33040lthioredoxin and At4g33050l 
embryo sac defective39; and At5gl3180INAC083 and 
At5gl3l90ILITAF domain. All these gene pairs were coor- 
dinate^ regulated, supporting the generalized observation 
that hlqltop6b can impact chromatin remodelling, leading 
to transcription of unrelated adjacent genes. As previously 
noted (Simkova et al, 2012) in the context of top6alcaa39 
characterization, there is in vivo molecular evidence that 
human TopoII(3 can generate transient double-stranded 
DNA breaks that permit exchanges in nucleosome-specific 
histone Hl-HMGB proteins, which promotes local changes 
of chromatin structure and leads to the transcriptional 
activation of target genes (Ju et al, 2006). Topoisomerase- 
mediated chromatin factor exchange represents an attractive 
mechanism for the transcriptional activation and repression 
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of different sets of genes. More broadly, the present findings 
also raise questions about the forces exerted by chromatin 
remodelling factors in shaping plant genomes during evo- 
lution, for example if TopoVI was acquired by lateral gene 
transfer from archaebacteria instead of having been lost from 
the animal and fungal kingdoms. The finding of statistically 
significant misregulation of histone genes in the hlq mutant 
(Supplementary Fig. SI 3, Datafile SI, Table S4 at JXB online) 
including Atlg06760l HISTONE HI is consistent with a chro- 
matin factor exchange model and implies that there may be a 
feedback mechanism for TopoVI to regulate partner histone 
expression. The present systems approach of transcriptome 
meta-analysis revealed that hlq mutants misregulate a subset 
of ontologically unrelated genes based on proximity, which 
may help to explain the pleiotropic phenotypes of Un3lhyp6l 
rhl3ltop6blhlq. 

Remarkably, mutants and inhibitors affecting the cell wall 
phenocopy several of the traits displayed in hlq. Plant cell 
walls are complex carbohydrate polymers, and it is estimated 
that ~10% of plant genomes are devoted to cell wall biogen- 
esis. However, knowledge of their biogenesis and functions 
in signal transduction of environmental stimuli and plant 
development is limited, even in model organisms. Starch is 
the major non-structural carbohydrate in plants. It serves as 
an important store of carbon that fuels plant metabolism and 
growth when they are unable to photosynthesize. Arabidopsis 
has proven to be a powerful genetic system for discovering 
how starch is synthesized and degraded (Streb and Zeeman, 
2012), and the present results (Figs 2, 6; Supplementary Fig. 
S9 at JXB online) suggest that HLQITOP6B can be added 
to the list of genes important for regulating starch and cell 
wall biosynthesis. However, caution should be exercised when 
inferring regulatory mechanisms because studies at the tran- 
scriptome level do not address post-translational regulation 
of metabolic pathways. Changes in sugar levels or other met- 
abolic intermediates might serve as the trigger for the fluxes 
into and out of the starch pool. It was proposed that macro- 
autophagy may play a role in degrading non-functional chlo- 
roplasts in mutants of starch metabolism (Stettler et al. , 2009). 
Hartung et al. (2002) also observed, at the ultrastructural 
level in the cytoplasm of top6blspoll-2 and top6alatspoll- 
3 mutants, agglomerations of protein bodies, peroxisomes, 
and small vacuoles indicative of autophagy. There is genetic 
evidence that accumulation of starch degradation intermedi- 
ates, such as disaccharides, causes severe cellular phenotypes 
such as chlorosis and chloroplast lysis, because mutations in 
fi-amylase2 and /3-amylase3, which are strongly down-regu- 
lated in hlq, prevent the release of maltose from the starch 
granule and decrease the severe cellular chlorosis observed 
in the maltose exporter 1 1 dpel double mutant (Stettler et al., 
2009). Upon export to the cytosol, maltose is metabolized 
by a glucosyl transfer reaction catalysed by the cytosolic dis- 
proportionating enzyme At2g40840/DPE2. Mutants of dpel 
accumulate starch and 200 times more maltose than the wild 
type, which causes severe chlorosis (Chia et al. , 2004; Stettler 
et al, 2009), similar to the observed hlq gene expression phe- 
notypes including significantly reduced DPE2 expression 
(f=0.009; Supplementary Datafile 1), and elevated expression 



of several autophagy genes (Supplementary Fig. Sll). The 
observation of significantly elevated expression of the two 
Glc6P/Pi transporter genes (At5g54800/GPT1 and Atlg61800l 
GPT2, f<0.003; Supplementary Datafile 1) is consistent with 
a model of starch metabolism intermediates functioning as 
regulators of carbon partitioning. Transcriptome profiling 
of Arabidopsis sucrose-inducible gene expression after 30 min 
relief from starvation also revealed similar affected processes 
to those reported here: up-regulation of GPT2 and differen- 
tial expression of 21 transcription regulators (43% concord- 
ance), 15 ubiquitin-targeting proteins (60% concordance, but 
opposite FC), trehalose phosphate synthases (50% concord- 
ance, but opposite FC), autophagy protein 8e (opposite FC), 
several glutaredoxins (75% concordance, plus one Atlg28480l 
GRX480, 7'<0.001), and carbon-scavenging enzymes for pro- 
line and myoinositol (all opposite FC: At4g39800IInositol- 
3-phosphate synthase! IPS1; At3g30775/proline oxidase; and 
myoinositol oxygenases MIOX2l4IAt2gl9800IAt4g26260) 
(Osuna et al. , 2007). That study also reported strong repres- 
sion of At2gl7880IDnaJ and of lipid catabolism genes after 
sucrose addition that mirrored the decrease of fatty acids in 
carbon-starved seedlings and their gradual recovery after 
sucrose addition. Of the handful of genes in Table 3 identi- 
fied from meta-analysis as highly significant, one is a down- 
regulated lipid biosynthetic gene and another is At2gl7880l 
DnaJ. Furthermore, the pervasive down-regulation of over- 
represented GDSL-motif lipases (P=0.01; Supplementary 
Fig. S12; Supplementary Datafile 1) suggests that lipids are 
being oxidized to fuel gluconeogenesis. These correlations 
strongly support the interpretation that hlq seedlings are 
starving, especially when taken together with results (Fig. 6, 
Table 2) that show that hlq leaves accumulate starch in the 
dark. Alternatively, the observed strong statistical association 
of hlq differentially expressed genes with protein degrada- 
tion (Supplementary Fig. S10) is consistent with the pheno- 
types of the detlcoplfus class of mutants disrupting the COP9 
signalosome, which is structurally and functionally similar to 
the 19S regulatory particle of the 26S proteosome and inter- 
acts with SCF-type E3 ubiquitin ligases to regulate targeted 
protein degradation (Smalle and Vierstra, 2004; Gusmaroli 
et al, 2007). 

AGP2, the regulatory subunit of the rate-limiting starch 
biosynthetic enzyme, is inactivated under oxidizing conditions 
(Streb and Zeeman, 2012), a state that the present transcrip- 
tome data and previous results (Simkova et al. ,2012) indicate 
exists in hlqltop6b and top6alcaa39 mutants. Soluble starch 
synthase 1 (At5g24300/SSl) has also been proposed to be 
regulated by the redox state (Glaring et al, 2012), and it was 
observed here to be up- regulated CP=0.005; Supplementary 
Datafile 1 at JXB online) in hlq, which is an incongruity com- 
pared with other starch biosynthesis genes. In potato, high 
levels of glucose and sucrose increase both the AGP activation 
state and starch synthesis. The signalling processes for AGP 
regulation have been proposed to occur through hexokinase 
for glucose and through SnRKl kinases for sucrose (Tiessen 
et al, 2003). HEXOKINA SE1 1 GL UC OSE-INSENSITIVE2 
was slightly up-regulated in hlq CP<0.05; Supplementary 
Datafile SI), whereas SnRKl kinases were not affected. 
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However, several of the 35 SnRK2/3 genes were signifi- 
cantly up-regulated (P<0.04; SnRK2.3; SnRK3.14/CIPK6/ 
SOS3-INTERACTING3; SnRK3.10/CIPK7; SnRK3.15/ 
CIPK14/PSK24; Supplementary Datafile SI) or down-regu- 
lated (P<0.008; SnRK2.1/ASK2; SnRK3.12/CIPK9/PKS6; 
SnRK3.17/CIPK3). Interestingly, these particular SnRKs 
have been shown to function in ABA, osmotic, and salt stress, 
flavonoid/sugar and calcium signalling, and auxin responses 
(Luan et al, 2002; Kolukisaoglu et al, 2004; Luan, 2009; 
Tripathi et al, 2009; Heyndrickx and Vandepoele, 2012), all 
processes described herein as affected by the hlq mutation 
and manifest in the pleiotropic phenotypes. Likewise, the hlq 
phenotypes of abnormal root hair development, sterility, and 
biotic stress responses through LRR-NBS receptor kinases is 
known to be controlled by calcium and ROS (Monshausen 
et al, 2008, 2009), processes significantly altered in hlq 
(Supplementary Fig. SI 1; Supplementary Datafile SI). 
Furthermore, there is evidence that the disaccharide treha- 
lose-6-phosphate (Tre6P) acts as a signalling intermediate in 
the sucrose -dependent activation of AGP (Lunn et al, 2006), 
which is strongly down-regulated in hlq (Supplementary 
Datafile SI). The finding that the trehalose marker PAD4 
and Tre6P biosynthetic gene TREHALOSE PHOSPHATE 
SYNTHASE11 were strongly up-regulated in hlq supports 
this model. However, the details of this mechanism remain to 
be elucidated, and the unexpected down-regulated expression 
of the AGP2 gene in hlq is not understood. Further analy- 
ses of trehalose and Tre6P in hlq and hypomorphic alleles 
could reveal the connections between the pathways of starch 
metabolism and the signalling pathways that regulate them, 
and provide insights into cell wall integrity-related mecha- 
nisms underlying environmental sensing, cell expansion, and 
morphogenesis (Hematy et al, 2007; Boisson-Dernier et al, 
201 1; Cheung and Wu, 201 1). 

The observed 1.9-fold up-regulation (P<0.0003) of 
At2g03770/ST1-Like sulfotransf erase, involved in brassinos- 
teroid metabolism (Marsolais et al. , 2007), is noted in connec- 
tion with bin3lhlq skotomorphogenesis phenotypes. Similarly, 
the hlq mutant showed significant down-regulation (f<0.02) 
of a brassinosteroid biosynthetic gene At5g05690lconstitu- 
tive photomorphogenic dwarflcabbage3lcyp90al. Mutants of 
cpdlcbb3lcyp90al display de-etiolation and de-repression of 
light-induced genes in the dark, dwarfism, male sterility, and 
activation of stress-regulated genes in the light. Likewise, the 
observed strong down-regulation in hlq (Table 5) of ABA1, 
CLA1, CED1, and NCED1 involved in ABA biosynthesis, 
skotomorphogenesis via cartenoid metabolism, and regulation 
of cell elongation and morphogenesis (Barrero et al. , 2008) is 
intriguing. Like hlq, mutants of cedllbodyguardl alpha-beta 
hydrolase are dwarfed, extremely sensitive to osmotic stress, 
and have abnormal leaves, collapsed cells, and reduced num- 
bers of trichomes (Wang et al, 2011). Taken together with 
the observed misregulation of lipid and isoprenoid metabo- 
lism pathways in hlq (Supplementary Fig. S9, Supplementary 
Datafile SI), it is speculated that TopoVI regulation of genes 
controlling wax or lipid metabolism may be the molecular 
mechanism underlying the unexpected genetic link (Wang 
et al. , 201 1 ) between cuticle biogenesis and ABA and osmotic 



stress responses. It is also speculated that TopoVI regulation 
of ATI GO 1 550/ B YPA SSI and At4g01360/BYPASS3 (signifi- 
cantly up-regulated in hlq; P<Q.Q0l, Supplementary Datafile 
SI) underlies the pleiotropic phenotypes of bpsl/3 mutants 
affecting root production of a mobile shoot morphogen of 
unknown structure derived from carotenoids (Lee et al, 
2012). 

The findings that hlqltop6b mutants resemble the detl 
coplfus class of mutants which accumulate the second- 
ary metabolites anthocyanins parallels several recent lines 
of genetic evidence that responses to environmental cues 
are tightly coupled through sensors that regulate carbon 
partitioning between primary (aromatic amino acid) and 
secondary (polyphenolic) metabolites (Tzin and Galili, 
2010). In addition to being a primary metabolite for pro- 
tein synthesis, phenylalanine availability is rate limiting for 
production of phenylpropanoids important for physiol- 
ogy, UV protectants, fragrances/flavours related to repro- 
duction, defence, hormonal pathways, and the synthesis 
of cell wall material that can comprise up to 50% of cap- 
tured photosynthetic carbon, depending on the species. 
One possible insight into the cell elongation, lignin and 
anthocyanin accumulation phenotypes of hlq (Fig. 3A, 
B; Supplementary Figs SIC, S5D at JXB online) is that 
stem elongation, biomass, and lignin biosynthesis in vas- 
cular bundles and fibre cells are controlled by stem- and 
root-specific expression of chloroplastic arogenate dehy- 
dratases At3g44720IADT4 and At5g22630IADT5 (Corea 
et al, 2012), which were significantly up-regulated in hlq 
mutants (7 > <0.0007; Supplementary Datafile 1). Another 
interesting facet is the link between ethylene and auxin 
signalling and transport, root hair development, cell elon- 
gation, and flavonols (Lewis et al, 2011): At2g20610l 
SUPERROOT1/HOOKLESS3/ROOTY is an auxin over- 
expressor encoding an aminotransferase that functions in 
cross-talk between auxin and ethylene and is a marker for 
glucosinolate signalling (Clay et al, 2009), and it was sig- 
nificantly down-regulated in hlq (P<0.01; Supplementary 
Datafile 1). Remarkably dc3 overexpressor2 {dorl) was iso- 
lated and shown that it was a new allele of surllrty from the 
same screen that produced hlq (Sun, 2003). Further analy- 
sis and classification of up- and down-regulated genes in 
other TopoVI subunit mutants and characterization at the 
genetic, protein, and molecular levels of a series of hypo- 
morphic TopoVI subunit B alleles (Fig. 5) may reveal the 
core function(s) of TopoVI. Studies on regulation of phe- 
nylalanine and polyphenolic biosynthetic pathways in both 
the chloroplast and cytosol, including plasma membrane- 
localized processes, should help unravel the complexities of 
lignin, cellulose, and polyphenolic contributions to struc- 
tural integrity and cellular functions at the transcriptional, 
metabolic, and tissue levels. 



Supplementary data 

Supplementary data are available at JXB online. 

Figure SI. Pleiotropic phenotypes of the hlq mutant. 
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Figure S2. Ectopic proDc3:GUS expression and ligni- 
fication in two cell- wall related mutants procustel-1 and 
boterol-1. 

Figure S3. Cell wall inhibitors tunicamycin and DCB phe- 
nocopy epidermal traits of hlq. 

Figure S4. Sanger sequencing results on both strands of 
independent amplicons spanning the hlq point mutation at 
chr3: 7267232, confirming the Illumina whole -genome re- 
sequencing result of G— >A transition. 

Figure S5. Dwarf, chlorotic, and anthocyanin-accumula- 
tion phenotypes of homozygotes for T-DNA insertion lines 
in At3g20780 and that disrupt exons 4 and 12. 

Figure S6. Pleiotropic root hair defects of anisotropic 
growth, bulging, and branching for T-DNA insertion lines in 
At3g20780/TOP6B. 

Figure S7. Non-complementation (segregating mutants) in 
Fj progeny of crosses between heterozygous hlql+ and het- 
erozygous T-DNA insertion lines, validated by PCR. 

Figure S8. Volcano plot of the hlq transcriptome profiling 
experiment, showing the relationship between fold change 
effects of the hlq genotype compared with the wild type and 
statistical significance based on three replicates for ~25 000 
genes on the microarray. 

Figure S9. Transcriptome profiling of hlq seedlings for 
general metabolic pathway effects using MAPMAN (Usadel 
et al, 2005). 

Figure S10. Transcriptome profiling of hlq seedlings for 
ubiquitin- and autophagy- dependent protein degradation 
pathways using MAPMAN. 

Figure SI 1. Transcriptome profiling of hlq seedlings to 
'Gene Regulation Overview' using MAPMAN. 

Figure SI 2. Transcriptome profiling of hlq seedlings to 
(A) 'Cellular Response Overview' and (B) 'Large Enzyme 
Families' pathways using MAPMAN. 

Figure SI 3. Transcriptome profiling of hlq seedlings to spe- 
cific families of transcription factors using MAPMAN. 

Datafile SI. Transcriptome Microarray log 2 FC of hlq/WT, 
with statistical significance for differential expression (sheet 
1). Sheet 2: overview of statistically significant over-repre- 
sented processes (Benjamini-Hochberg corrected) affected by 
hlq, based on Wilcoxon Rank Sum Test of the transcriptome 
data set calculated by MAPMAN. Sheet 3: list of genes in 
statistically over-represented bins. 

Data file S2. Detailed characterization of the hlq 
transcriptome. 

Table SI. Cell length parameters of wild-type and hlq 
mutant hypocotyl and roots. 

Table S2 (Datafile 1, sheet 4). Analysis of the 1079 most 
differentially expressed genes in the hlq mutant for chro- 
mosomal adjacency and co-regulation, compared with a 
comparable number of 'control' genes with the smallest FC 
symmetrically distributed about FC=0. 

Table S3 (Datafile 1, sheet 5). Concordance between hlq 
transcriptome results and previously published lists of mis- 
regulated genes in bin3ltop6b (Yin et al, 2002) and top6al 
caa39 (Simkova et al, 2012). 

Table S4. List of HISTONE, WRKY, and AUX/IAA genes 
significantly mis-regulated in hlq mutant. 



Table S5 (Datafile 1, sheet 6). Analysis of concord- 
ance between hlq misregulated genes compared with 
published transcriptome data sets for ABA RESPONSI 
VE17-(PYRABACTIN-RESISTANCE-LIKE12 homolo- 
gous) overexpressing transgenic plants (Krishnaswamy 
et al, 2008) and ABA-regulated genes (Matsui et al, 
2008). 

Table S6. Primers for new markers developed during fine- 
mapping the hlq mutant, and other primers. 

Acknowledgements 

The authors thank Terry Thomas for providing the proDc3:GUS reporter 
lines 7-2 and 11-1, Ralph Quatrano for support, Lucy Smith for back- 
crossing hlq, the SEM facility at the Department of Pathology, University 
of Hong Kong for technical assistance with scanning electron micros- 
copy, Thomas Altmann, Shibo Li, and Megan Sweeney for assistance 
with microarrays, Bruce Roe and Scott Dowd for preliminary sequencing 
and data analysis, respectively, Ruobai Sun for assistance with Illumina 
short read data sets, an anonymous reviewer for helpful suggestions on 
the manuscript, and the TTU High Performance Computing Center for 
support of the Hrothgar computer cluster. This work was supported by 
grant #2006-34186-16976-sub-J55-Q01097 from the USDA Southwest 
Consortium for Plant Genetics and Water Resources to CDR and DWG. 
The funders had no role in the study design, in the collection, analysis, or 
interpretation of data, or in the writing of the manuscript or decision to 
submit the manuscript for publication. 



References 

Alonso-Simon A, Garcia-Angulo P, Melida H, Encina A, Alvarez JM, 
Acebes JL. 201 1 . The use of FTIR spectroscopy to monitor modifications 
in plant cell wall architecture caused by cellulose biosynthesis inhibitors. 
Plant Signaling and Behavior 6, 1 1 04-1 1 1 0. 

Arioli T, Peng L, Betzner AS, ef al. 1 998. Molecular analysis of cellulose 
biosynthesis in Arabidopsis. Science 279, 717-720. 

Balasubramanian R. 2003. Harlequin (hlq): an Arabidopsis mutant 
that ectopically expresses Dc3-GUS and shows defects in cell wall 
morphogenesis. PhD Thesis, Hong Kong University of Science and 
Technology, http://ustlib.ust.hk/. 

Barrero JM, Rodriguez PL, Quesada V, Alabadi D, Blazquez MA, 
Boutin JP, Marion-Poll A, Ponce MR, Micol JL. 2008. TbeABM 
gene and carotenoid biosynthesis are required for late skotomorphogenic 
growth in Arabidopsis thaliana. Plant, Cell and Environment 31, 227-234. 

Bell CJ, Ecker JR. 1 994. Assignment of 30 microsatellite loci to the 
linkage map of Arabidopsis. Genomics 19, 137-144. 

Bichet A, Desnos T, Turner S, Grandjean O, Hofte H. 2001 . BOTEROI 
is required for normal orientation of cortical microtubules and anisotropic 
cell expansion in Arabidopsis. The Plant Journal 25, 1 37-1 48. 

Boisson-Dernier A, Kessler SA, Grossniklaus U. 201 1 . The walls have 
ears: the role of plant CrRLKI Ls in sensing and transducing extracellular 
signals. Journal of Experimental Botany 62, 1 581 -1 591 . 

Bougourd S, Marrison J, Haseloff J. 2000. An aniline blue staining 
procedure for confocal microscopy and 3D imaging of normal and 
perturbed cellular phenotypes in mature Arabidopsis embryos. The Plant 
Journal 24, 543-550. 

Breuer C, Stacey NJ, West CE, Zhao Y, Chory J, Tsukaya H, Azumi 
Y, Maxwell A, Roberts K, Sugimoto-Shirasu K. 2007. BIN4, a novel 
component of the plant DNA topoisomerase VI complex, is required for 
endoreduplication in Arabidopsis. The Plant Cell 19, 3655-3668. 

Brocard-Gifford I, Lynch TJ, Garcia ME, Malhotra B, Finkelstein 

RR. 2004. The Arabidopsis thaliana ABSCISIC ACID-INSENSITIVE8 locus 
encodes a novel protein mediating abscisic acid and sugar responses 
essential for growth. The Plant Cell 16, 406-421 . 

Brou C, Kuhn A, Staub A, Chaudhary S, Grummt I, Davidson I, Tora 

L. 1993. Sequence-specific transactivators counteract topoisomerase 
ll-mediated inhibition of in vitro transcription by RNA polymerases I and II. 
Nucleic Acids Research 21 , 401 1 -401 8. 



TopoVI regulation of growth and development | 4237 



Brownleader MD, Hopkins J, Mobasheri A, Dey PM, Jackson P, 
Trevan M. 2000. Role of extensin peroxidase in tomato (Lycopersicon 
esculentum Mill.) seedling growth. Planta 210, 668-676. 

Burk DH, Liu B, Zhong R, Morrison WH, Ye Z-H. 2001 . A katanin-like 
protein regulates normal cell wall biosynthesis and cell elongation. The 
Plant Cell A3, 807-827. 

Cafio-Delgado Al, Metzlaff K, Bevan MW. 2000. The eli1 mutation 
reveals a link between cell expansion and secondary cell wall formation in 
Arabidopsis thaliana. Development 127, 3395-3405. 

Chak RKF, Thomas TL, Quatrano RS, Rock CD. 2000. The genes 
ABU and ABI2 are involved in abscisic acid- and drought-inducible 
expression of the Daucus carota Dc3 promoter in guard cells of transgenic 
Arabidopsis thaliana (L.) Heynh. Planta 210, 875-883. 

Chen X-Y, Kim J-Y. 2009. Callose synthesis in higher plants. Plant 
Signaling and Behavior 4, 489-492. 

Chen ZZ, Hong XH, Zhang HR, Wang YQ, Li X, Zhu JK, Gong ZZ. 

2005. Disruption of the cellulose synthase gene, AtCesA8/IRX1 , enhances 
drought and osmotic stress tolerance in Arabidopsis. The Plant Journal 43, 
273-283. 

Cheng J-C, Lertpiriyapong K, Wang S, Sung ZR. 2000. The role of the 
Arabidopsis ELD1 gene in cell development and photomorphogenesis in 
darkness. Plant Physiology 123, 509-520. 

Cheung AY, Wu H-M. 201 1 . THESEUS 1 , FERONIA and relatives: a 
family of cell wall-sensing receptor kinases? Current Opinion in Plant 
Biology 1 4, 632-641. 

Chia T, Thorneycroft D, Chappie A, Messerli G, Chen J, Zeeman SC, 
Smith SM, Smith AM. 2004. A cytosolic glucosyltransferase is required 
for conversion of starch to sucrose in Arabidopsis leaves at night. The 
Plant Journal 37, 853-863. 

Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land 
SJ, Lu X, Ruden DM. 2012. A program for annotating and predicting 
the effects of single nucleotide polymorphisms, SnpEff: SNPs in the 
genome of Drosophila melanogaster strain win 8; iso-2; iso-3. Fly 6, 
80-92. 

Clay NK, Adio AM, Denoux C, Jander G, Ausubel FM. 2009. 
Glucosinolate metabolites required for an Arabidopsis innate immune 
response. Science 323, 95-1 01 . 

Corbett KD, Berger JM. 2006. Structural basis for topoisomerase VI 
inhibition by the anti-Hsp90 drug radicicol. Nucleic Acids Research 34, 
4269-4277. 

Corea ORA, Ki C, Cardenas CL, Kim S-J, Brewer SE, Patten 
AM, Davin LB, Lewis NG. 2012. Arogenate dehydratase isoenzymes 
profoundly and differentially modulate carbon flux into lignins. Journal of 
Biological Chemistry 287, 1 1 446-1 1 459. 

Desikan R, Hagenbeek D, Neill SJ, Rock CD. 1 999. Flow cytometry 
and surface plasmon resonance analyses demonstrate that the monoclonal 
antibody JIM19 interacts with a rice cell surface component involved in 
abscisic acid signalling in protoplasts. FEBS Letters 456, 257-262. 

Ding L, Zhu J-K. 1997. A role for arabinogalactan-proteins in root 
epidermal cell expansion. Planta 203, 289-294. 

Evans-Roberts KM, Breuer C, Wall MK, Sugimoto-Shirasu K, 
Maxwell A. 2010. Arabidopsis thaliana GYRB3 does not encode a DNA 
gyrase subunit. PLoS One 5, e9899. 

Fagard M, Desnos T, Desprez T, Goubet F, Refregier G, Mouille G, 
McCann M, Rayon C, Vernhettes S, Hofte H. 2000. PROCUSTE1 
encodes a cellulose synthase required for normal cell elongation 
specifically in roots and dark-grown hypocotyls of Arabidopsis. The Plant 
Ce//12, 2409-2423. 

Finkelstein RR, Gampala SSL, Rock CD. 2002. Abscisic acid signaling 
in seeds and seedlings. The Plant Cell 14, S15-S45. 

Fry SC. 2000. The growing plant cell wall: chemical and metabolic 
analysis , reprint edn. Caldwell, NJ: The Blackburn Press. 

Gadjev I, Vanderauwera S, Gechev TS, Laloi C, Minkov IN, Shulaev 
V, Apel K, Inze D, Mittler R, Van Breusegem F. 2006. Transcriptomic 
footprints disclose specificity of reactive oxygen species signaling in 
Arabidopsis. Plant Physiology 141, 436-445. 

Galbraith DW, Harkins KR, Knapp S. 1 991 . Systemic endopolyploidy in 
Arabidopsis thaliana. Plant Physiology 96, 985-989. 

Galbraith DW, Janda J, Lambert GM. 201 1 . Multiparametric 
analysis, sorting, and transcriptional profiling of plant protoplasts 



and nuclei according to cell type. Methods in Molecular Biology 699, 
407-429. 

Gao M, Showalter AM. 1 999. Yariv reagent treatment induces 
programmed cell death in Arabidopsis cell cultures and implicates 
arabinogalactan protein involvement. The Plant Journal 19, 321-331 . 

Gentleman R, Carey V, Bates D, etal. 2004. Bioconductor: open 
software development for computational biology and bioinformatics. 
Genome Biology 5, R80. 

Glaring MA, Skryhan K, Kotting O, Zeeman SC, Blennow A. 2012. 
Comprehensive survey of redox sensitive starch metabolising enzymes in 
Arabidopsis thaliana. Plant Physiology and Biochemistry 58, 89-97. 

Guan Y, Nothnagel EA. 2004. Binding of arabinogalactan proteins by 
Yariv phenylglycoside triggers wound-like responses in Arabidopsis cell 
cultures. Plant Physiology 135, 1346-1366. 

Gusmaroli G, Figueroa P, Serino G, Deng XW. 2007. Role of the MPN 
subunits in COP9 signalosome assembly and activity, and their regulatory 
interaction with Arabidopsis Cullin3-based E3 ligases. The Plant Cell 19, 
564-581 . 

Hematy K, Sado P-E, Van Tuinen A, Rochange S, Desnos T, 
Balzergue S, Pelletier S, Renou J-P, Hofte H. 2007. A receptor-like 
kinase mediates the response of Arabidopsis cells to the inhibition of 
cellulose synthesis. Current Biology 17, 922-931 . 

Hartung F, Angelis KJ, Meister A, Schubert I, Melzer M, Puchta H. 2002. 
An archaebacterial topoisomerase homolog not present in other eukaryotes is 
indispensable for cell proliferation of plants. Current Biology 12, 1 787-1 791 . 

Hartung F, Puchta H. 2001 . Molecular characterization of homologues of 
both subunits A (SP01 1) and B of the archaebacterial topoisomerase 6 in 
plants. Gene 271, 81-86. 

Heyndrickx KS, Vandepoele K. 2012. Systematic identification of 
functional plant modules through the integration of complementary data 
sources. Plant Physiology 159, 884-901 . 

His I, Driouich A, Nicol F, Jauneau A, Hofte H. 2001 . Altered pectin 
composition in primary cell walls of korrigan, a dwarf mutant of Arabidopsis 
deficient in a membrane-bound endo-1 ,4-p-glucanase. Planta 212, 
348-358. 

Hu W, Ma H. 2006. Characterization of a novel putative zinc finger 
gene MIF1: involvement in multiple hormonal regulation of Arabidopsis 
development. The Plant Journal 45, 399-422. 

Huh SU, Lee SB, Kim HH, Paek KH. 2012. ATAF2, a NAC transcription 
factor, binds to the promoter and regulates NIT2 gene expression involved 
in auxin biosynthesis. Molecules and Cells 34, 305-313. 

Jain M, Tyagi AK, Khurana JP. 2006. Overexpression of putative 
topoisomerase 6 genes from rice confers stress tolerance in transgenic 
Arabidopsis plants. FEBS Journal 273, 5245-5260. 

Jain M, Tyagi AK, Khurana JP. 2008. Constitutive expression of a 
meiotic recombination protein gene homolog, OsTOP6A1 , from rice 
confers abiotic stress tolerance in transgenic Arabidopsis plants. Plant Cell 
Reports 27, 767-778. 

Jakoby MJ, Falkenhan D, Mader MT, Brininstool G, Wischnitzki 
E, Platz N, Hudson A, Hulskamp M, Larkin J, Schnittger A. 2008. 
Transcriptional profiling of mature Arabidopsis trichomes reveals that 
NOECK encodes the MIXTA-Like transcriptional regulator MYB106. Plant 
Physiology 148, 1583-1602. 

Jander G, Norris SR, Rounsley SD, Bush DF, Levin IM, Last RL. 

2002. Arabidopsis map-based cloning in the post-genome era. Plant 
Physiology 129, 440-450. 

Jiang T, Zhang X-F, Wang X-F, Zhang D-P. 201 1 . Arabidopsis 
3-Ketoacyl-CoA Thiolase-2 (KAT2), an enzyme of fatty acid beta-oxidation, 
is involved in ABA signal transduction. Plant and Cell Physiology 52, 
528-538. 

Ju B-G, Lunyak W, Perissi V, Garcia-Bassets I, Rose DW, Glass 
CK, Rosenfeld MG. 2006. A topoisomerase lib-mediated dsDNA break 
required for regulated transcription. Science 312, 1798-1802. 

Kirik V, Schrader A, Uhrig JF, Hulskamp M. 2007. MIDGET unravels 
functions of the Arabidopsis Topoisomerase VI complex in DNA 
endoreduplication, chromatin condensation, and transcriptional silencing. 
The Plant Cell 19, 31 00-31 1 0. 

Kohler A, Schwindling S, Conrath U. 2000. Extraction and quantitative 
determination of callose from Arabidopsis leaves. Biotechniques 28, 
1084-1086. 



4238 I Mittal ef a/. 



Kolukisaoglu U, Weinl S, Blazevic D, Batistic O, Kudla J. 2004. 
Calcium sensors and their interacting protein kinases: genomics of the 
Arabidopsis and rice CBL-CIPK signaling networks. Plant Physiology 134, 
43-58. 

Komalavilas P, Zhu JK, Nothnagel EA. 1 991 . Arabinogalactan-proteins 
from the suspension-culture medium and plasma-membrane of rose cells. 
Journal of Biological Chemistry 266, 1 5956-1 5965. 

Kong D, Karve R, Willet A, Chen M-K, Oden J, Shpak ED. 2012. 
Regulation of plasmodesmatal permeability and stomatal patterning by the 
glycosyltransferase-like protein KOBITOI . Plant Physiology 1 59, 1 56-1 68. 

Konieczny A, Ausubel FM. 1 993. A procedure for mapping Arabidopsis 
mutations using co-dominant ecotype-specific PCR-based markers. The 
Plant Journal 4, 403-41 0. 

Kosak ST, Scalzo D, Alworth SV, Li F, Palmer S, Enver T, Lee JSJ, 
Groudine M. 2007. Coordinate gene regulation during hematopoiesis is 
related to genomic organization. PLoS Biology 5, e309. 

Kurata T, Yamamoto KT. 1998. petitl , a conditional growth mutant of 
Arabidopsis defective in sucrose-dependent elongation growth. Plant 
Physiology 118, 793-801 . 

Lee D-K, Van Norman JM, Murphy C, Adhikari E, Reed JW, Sieburth 

LE. 201 2. In the absence of BYPASS1 -related gene function, the bps 
signal disrupts embryogenesis by an auxin-independent mechanism. 
Development 139, 805-815. 

Lee MM, Schiefelbein J. 1 999. WEREWOLF, a MYB-related protein in 
Arabidopsis, is a position-dependent regulator of epidermal cell patterning. 
Cell 99, 473-483. 

Lewis DR, Ramirez MV, Miller ND, Vallabhaneni P, Ray WK, Helm 
RF, Winkel BSJ, Muday GK. 201 1 . Auxin and ethylene induce flavonol 
accumulation through distinct transcriptional networks. Plant Physiology 
156, 144-164. 

Li H, Durbin R. 2010. Fast and accurate long-read alignment with 
Burrows-Wheeler transform. Bioinformatics 26, 589-595. 

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth 
G, Abecasis G, Durbin R, 1000 Genome Project Data Processing 
Subgroup. 2009. The sequence alignment/map format and SAMtools. 
Bioinformatics 25, 2078-2079. 

Liu Z-Q, Yan L, Wu Z, Mei C, Lu K, Yu Y-T, Liang S, Zhang X-F, Wang 
X-F, Zhang D-P. 201 2. Cooperation of three WRKY-domain transcription 
factors WRKY1 8, WRKY40, and WRKY60 in repressing two ABA- 
responsive genes ABI4 and ABI5 in Arabidopsis. Journal of Experimental 
Botany 63, 6371-6392. 

Lu B, Sun W, Zhang S, Zhang C, Qian J, Wang X, Gao R, Dong H. 

201 1 . HrpN Ea -induced deterrent effect on phloem feeding of the green 
peach aphid Myzus persicae requires AtGSL5 and AtMYB44 genes in 
Arabidopsis thaliana. Journal of Bioscience 36, 123-137. 

Luan S. 2009. The CBL-CIPK network in plant calcium signaling. Trends 
in Plant Science 14, 37-42. 

Luan S, Kudla J, Rodriguez-Concepcion M, Yalovsky S, Gruissem 

W. 2002. Calmodulins and calcineurin-like proteins: calcium sensors for 
specific signal response coupling in plants. The Plant Cell 14, S389-S400. 

Lukowitz W, Nickle TC, Meinke DW, Last RL, Conklin PL, 
Somerville CR. 2001 . Arabidopsis cyt 1 mutants are deficient in a 
mannose-1 -phosphate guanylyltransferase and point to a requirement 
of N-linked glycosylation for cellulose biosynthesis. Proceedings of 
the National Academy of Sciences, USA 98, 2262-2267. 

Lunn JE, Feil R, Hendriks JHM, Gibon Y, Morcuende R, Osuna D, 
Scheible W-R, Carillo P, Hajirezaei M-R, Stitt M. 2006. Sugar-induced 
increases in trehalose 6-phosphate are correlated with redox activation 
of ADPglucose pyrophosphorylase and higher rates of starch synthesis in 
Arabidopsis thaliana. Biochemical Journal 397, 139-148. 

Luo Q-J, Mittal A, Jia F, Rock C. 2012. An autoregulatory feedback 
loop involving PAP1 and TAS4 in response to sugars in Arabidopsis. Plant 
Molecular Biology 80, 1 1 7-1 29. 

Makarevitch I, Somers DA. 2006. Association of 
Arabidopsis topoisomerase 1 1 A cleavage sites with functional 
genomic elements and T-DNA loci. The Plant Journal 48, 
697-709. 

Marks MD, Wenger JP, Gilding E, Jilk R, Dixon RA. 2009. 
Transcriptome analysis of Arabidopsis wild-type and gl3-sst sim trichomes 



identifies four additional genes required for trichome development. 
Molecular Plant 2, 803-822. 

Marsolais F, Boyd J, Paredes Y, Schinas A-M, Garcia M, Elzein 
S, Varin L. 2007. Molecular and biochemical characterization of two 
brassinosteroid sulfotransferases from Arabidopsis, AtST4a (At2g14920) 
and AtST1 (At2g03760). Planta 225, 1233-1244. 

McCann MC, Hammouri M, Wilson R, Belton P, Roberts K. 1992. 
Fourier transform infrared microspectroscopy is a new way to look at plant 
cell walls. Plant Physiology 100, 1940-1947. 

McNamara S, Wang H, Hanna N, Miller WH. 2008. Topoisomerase lib 
negatively modulates retinoic acid receptor a function: a novel mechanism 
of retinoic acid resistance. Molecular and Cellular Biology 28, 2066-2077. 

Mochizuki N, Brusslan JA, Larkin R, Nagatani A, Chory J. 2001 . 
Arabidopsis genomes uncoupled 5 (GUN5) mutant reveals the involvement 
of Mg-chelatase H subunit in plastid-to-nucleus signal transduction. 
Proceedings of the National Academy of Sciences, USA 98, 2053-2058. 

Monshausen GB, Bibikova TN, Weisenseel MH, Gilroy S. 2009. 
Ca 2+ regulates reactive oxygen species production and pH during 
mechanosensing in Arabidopsis roots. The Plant Cell 21, 2341-2356. 

Monshausen GB, Messerli MA, Gilroy S. 2008. Imaging of the Yellow 
Cameleon 3.6 indicator reveals that elevations in cytosolic Ca 2+ follow 
oscillating increases in growth in root hairs of Arabidopsis. Plant Physiology 
147, 1690-1698. 

Nemhauser JL, Hong FX, Chory J. 2006. Different plant hormones 
regulate similar processes through largely nonoverlapping transcriptional 
responses. Cell 126, 467-475. 

Osuna D, Usadel B, Morcuende R, Gibon Y, ef al. 2007. Temporal 
responses of transcripts, enzyme activities and metabolites after adding 
sucrose to carbon-deprived Arabidopsis seedlings. The Plant Journal 49, 
463-491 . 

Pagant S, Bichet A, Sugimoto K, Lerouxel O, Desprez T, McCann 
M, Lerouge P, Vernhettes S, Hofte H. 2002. KOBIT01 encodes a novel 
plasma membrane protein necessary for normal synthesis of cellulose 
during cell expansion in Arabidopsis. The Plant Cell 14, 2001-2013. 

Paquette J, Tokuyasu T. 2010. EGAN: exploratory gene association 
networks. Bioinformatics 26, 285-286. 

Park MH, Suzuki Y, Chono M, Knox JP, Yamaguchi I. 2003. CsAGPI , 
a gibberellin-responsive gene from cucumber hypocotyls, encodes a 
classical arabinogalactan protein and is involved in stem elongation. Plant 
Physiology 1 31 , 1 450-1 459. 

Rautengarten C, Ebert B, Herter T, Petzold CJ, Ishii T, 
Mukhopadhyay A, Usadel B, Scheller HV. 201 1 . The interconversion of 
UDP-arabinopyranose and UDP-arabinofuranose is indispensable for plant 
development in Arabidopsis. The Plant Cell 23, 1 373-1 390. 

Rock CD, Sakata Y, Quatrano RS. 2010. Stress signaling I: the role 
of abscisic acid. In: Pareek A, Sopory SK, Bohnert HJ, Govindjee, eds. 
Abiotic stress adaptation in plants: physiological, molecular and genomic 
foundation. Dordrecht, The Netherlands: Springer, 33-73. 

Rock CD, Sun X. 2005. Crosstalk between ABA and auxin signaling 
pathways in roots of Arabidopsis thaliana (L.) Heynh. Planta 222, 98-106. 

Schiefelbein JW, Somerville C. 1 990. Genetic control of root hair 
development in Arabidopsis thaliana. The Plant Cell 2, 235-243. 

Schneeberger K, Ossowski S, Ott F, ef al. 201 1 . Reference-guided 
assembly of four diverse Arabidopsis thaliana genomes. Proceedings of 
the National Academy of Sciences, USA 108, 10249-10254. 

Schneider K, Mathur J, Boudonck K, Wells B, Dolan L, Roberts K. 

1998. The ROOT HAIRLESS 1 gene encodes a nuclear protein required 
for root hair initiation in Arabidopsis. Genes and Development 12, 
2013-2021. 

Schrader A, Welter B, Hulskamp M, Hoecker U, Uhrig JF. 2013. 
MIDGET connects COP1 -dependent development with endoreduplication 
in Arabidopsis thaliana. The Plant Journal 75, 67-79. 

Selvendran RR, O'Neill MA. 1987. Isolation and analysis of cell walls 
from plant material. Methods in Biochemical Analysis 32, 25-1 53. 

Shang Y, Yan L, Liu Z-Q, ef al. 201 0. The Mg-chelatase H subunit of 
Arabidopsis antagonizes a group of WRKY transcription repressors to 
relieve ABA-responsive genes of inhibition. The Plant Cell 22, 1 909-1 935. 

Shen Y-Y, Wang X-F, Wu F-Q, ef al. 2006. The Mg-chelatase H subunit 
is an abscisic acid receptor. Nature 443, 823-826. 



TopoVI regulation of growth and development | 4239 



Shevell DE, Leu W-M, Gillmor CS, Xia G, Feldmann KA, Chua N-H. 

1994. EMB30 is essential for normal cell division, cell expansion, and cell 
adhesion in Arabidopsis and encodes a protein that has similarity to Sec7. 
Cell 77, 1051-1062. 

Simkova K, Moreau F, Pawlak P, Vriet C, Baruah A, Alexandre 
C, Hennig L, Apel K, Laloi C. 2012. Integration of stress-related 
and reactive oxygen species-mediated signals by Topoisomerase VI in 
Arabidopsis thaliana. Proceedings of the National Academy of Sciences, 
USA 109, 16360-16365. 

Smalle J, Vierstra RD. 2004. The ubiquitin 26S proteasome proteolytic 
pathway. Annual Review of Plant Biology 55, 555-590. 

Smyth GK, Michaud J, Scott HS. 2005. Use of within-array replicate 
spots for assessing differential expression in microarray experiments. 
Bioinformatics 21, 2067-2075. 

Stettler M, Eicke S, Mettler T, Messerli G, Hortensteiner S, Zeeman 

SC. 2009. Blocking the metabolism of starch breakdown products in 
Arabidopsis leaves triggers chloroplast degradation. Molecular Plant 2, 
1233-1246. 

Streb S, Zeeman SC. 201 2. Starch metabolism in Arabidopsis. The 
Arabidopsis Book. 10, e0160. 

Subramanian S, Rajagopal B, Rock CD. 2002. Harlequin (hlq) and 
short blue root (sbr), two Arabidopsis mutants that ectopically express an 
abscisic acid- and auxin-inducible transgenic carrot promoter and have 
pleiotropic effects on morphogenesis. Plant Molecular Biology 49, 93-105. 

Sugimoto-Shirasu K, Roberts GR, Stacey NJ, McCann MC, Maxwell 
A, Roberts K. 2005. RHL1 is an essential component of the plant DNA 
topoisomerase VI complex and is required for ploidy-dependent cell 
growth. Proceedings of the National Academy of Sciences, USA 102, 
18736-18741. 

Sugimoto-Shirasu K, Stacey NJ, Corsar J, Roberts K, McCann 

MC. 2002. DNA topoisomerase VI is essential for endoreduplication in 
Arabidopsis. Current Biology 12, 1 782-1 786. 

Sun X. 2003. Characterization of an auxin- and abscisic acid-inducible 
reporter gene: Dc3-GUS in reported auxin mutants, and mutant 
screening based on auxin responsive Dc3-GUS expression. MSc thesis. 
Hong Kong University of Science and Technology, http://ustlib.ust.hk/. 

Susek RE, Ausubel FM, Chory J. 1993. Signal transduction mutants 
of Arabidopsis uncouple nuclear CAB and RBCS gene expression from 
chloroplast development. Cell 74, 787-799. 

Szklarczyk D, Franceschini A, Kuhn M, ef a/. 201 1 . The STRING 
database in 201 1 : functional interaction networks of proteins, globally 
integrated and scored. Nucleic Acids Research 39, D561-D568. 

Tiessen A, Prescha K, Branscheid A, Palacios N, McKibbin R, 
Halford NG, Geigenberger P. 2003. Evidence that SNF1 -related 
kinase and hexokinase are involved in separate sugar-signalling 
pathways modulating post-translational redox activation of ADP-glucose 
pyrophosphorylase in potato tubers. The Plant Journal 35, 490-500. 

Ton J, Jakab G, Toquin V, Flors V, lavicoli A, Maeder MN, Metraux 
JP, Mauch-Mani B. 2005. Dissecting the beta-aminobutyric acid-induced 
priming phenomenon in Arabidopsis. The Plant Cell 17, 987-999. 

Ton J, Mauch-Mani B. 2004. beta-amino-butyric acid-induced resistance 
against necrotrophic pathogens is based on ABA-dependent priming for 
callose. The Plant Journal 38, 1 1 9-1 30. 



Tripathi V, Parasuraman B, Laxmi A, Chattopadhyay D. 2009. CIPK6, 
a CBL-interacting protein kinase is required for development and salt 
tolerance in plants. The Plant Journal 58, 778-790. 

Tsuzuki T, Takahashi K, Inoue S-l, Okigaki Y, Tomiyama M, Hossain 
M, Shimazaki K-l, Murata Y, Kinoshita T. 201 1 . Mg-chelatase H 
subunit affects ABA signaling in stomatal guard cells, but is not an ABA 
receptor in Arabidopsis thaliana. Journal of Plant Research 124, 527-538. 

Tzin V, Galili G. 2010. New insights into the shikimate and aromatic 
amino acids biosynthesis pathways in plants. Molecular Plant 3, 956-972. 

Usadel B, Nagel A, Steinhauser D, et al. 2006. PageMan: an interactive 
ontology tool to generate, display, and annotate overview graphs for 
profiling experiments. BMC Bioinformatics 7, 535. 

Usadel B, Nagel A, Thimm O, et al. 2005. Extension of the 
visualization tool MapMan to allow statistical analysis of arrays, display 
of corresponding genes, and comparison with known responses. Plant 
Physiology 138, 1 1 95-1 204. 

Varga-Weisz PD, Wilm M, Bonte E, Dumas K, Mann M, Becker PB. 

1 997. Chromatin-remodelling factor CHRAC contains the ATPases ISWI 
and topoisomerase II. Nature 388, 598-602. 

Vieira P, Kyndt T, Gheysen G, Engler JD. 2013. An insight into critical 
endocycle genes for plant-parasitic nematode feeding sites establishment. 
Plant Signaling and Behavior 8, e24223. 

Vos SM, Tretter EM, Schmidt BH, Berger JM. 201 1 . All tangled up: 
how cells direct, manage and exploit topoisomerase function. Nature 
Reviews Molecular Cell Biology 12, 827-841 . 

Wang Z-Y, Xiong L, Li W, Zhu J-K, Zhu J. 201 1 . The plant cuticle is 
required for osmotic stress regulation of abscisic acid biosynthesis and 
osmotic stress tolerance in Arabidopsis. The Plant Cell 23, 1971-1984. 

Willats WGT, Knox JP. 1996. A role for arabinogalactan-proteins in plant 
cell expansion: evidence from studies on the interaction of p-glucosyl 
Yariv reagent with seedlings of Arabidopsis thaliana. The Plant Journal 9, 
919-925. 

Wu F-Q, Xin Q, Cao Z, et al. 2009. The magnesium-chelatase H subunit 
binds abscisic acid and functions in abscisic acid signaling: new evidence 
in Arabidopsis. Plant Physiology 150, 1940-1954. 

Xie S, Lam E. 1994. Abundance of nuclear DNA topoisomerase II 
is correlated with proliferation in Arabidopsis thaliana. Nucleic Acids 
Research 22, 5729-5736. 

Xin Z, Chen J. 2012. A high throughput DNA extraction method with high 
yield and quality. Plant Methods 8, 26. 

Xu Y-H, Liu R, Yan L, Liu Z-Q, Jiang S-C, Shen Y-Y, Wang X-F, 
Zhang D-P. 2012. Light-harvesting chlorophyll a/b-binding proteins are 
required for stomatal response to abscisic acid in Arabidopsis. Journal of 
Experimental Botany 63, 1 095-1 1 06. 

Yin Y, Cheong H, Friedrichsen D, Zhao Y, Hu J, Mora-Garcia S, 
Chory J. 2002. A crucial role for the putative Arabidopsis topoisomerase 
VI in plant growth and development. Proceedings of the National Academy 
of Sciences, USA 99, 10191-10196. 

York WS, Darvill AG, McNeil M, Stevenson TT, Albersheim P. 1 986. 
Isolation and characterization of plant cell walls and cell wall components. 
Methods in Enzymology 118, 3-40. 

Zhong R, Ripperger A, Ye Z-H. 2000. Ectopic deposition of lignin in the 
pith of stems of two Arabidopsis mutants. Plant Physiology 123, 59-70. 



